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SUMMARY 


This report presents data from which one may obtain the probability that a 
pulsed-type radar system will detect a given target at any range. This is in con- 
trast to the usual method of obtaining radar range as a single number, which can be 
taken mathematically to imply that the probability of detection is zero at any range 
greater than this number, and certainty at any range less than this number. 


Three variables, which have so far received little quantitative attention in the 
subject of radar range, are introduced in the theory: 


1. The time taken to detect the target. 


2. The average time interval between false alarms 
(false indications of targets). 


3. The number of pulses integrated. 


It is shown briefly how the results for pulsed-type systems may be applied in the 
case of continuous-wave systems. 


Those concerned with systems analysis problems including radar performance may 
profitably use this work as one link in a chain involving several probabilities. For 
instance, the probability that a given aircraft will be detected at least once while 
flying any given path through a specified model radar network may be calculated using 
the data presented here as a basis, provided that additional probability data on such 
things as outage time etc., are available. 


The theory developed here does not take account of interference such as clutter 
or man-made static, but assumes only random noise at the receiver input. Also, an 
ideal type of electronic integrator and detector are assumed. Thus the results are 
the best that can be obtained under ideal conditions. It is not too difficult, how- 
ever, to make reasonable assumptions which will permit application of the results to 
the currently available types of radar. 


The first part of this report is a restatement of known radar fundamentals and 
supplies continuity and background for what follows. 


The mathematical part of the theory is not contained herein, but will be issued 
subsequently as a separate report ‘?*), 
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SYMBOLS 


A, = effective area of antenna for receiving. 

B = beamwidth of antenna. 

c = velocity of light. 

C = total shunt capacity of input circuit. 

) = factor which accounts for losses in transmission lines, T R switches, atmos- 


pheric absorption, etc. 
= rms noise voltage. 
= transmitted energy per pulse. 
= received energy per pulse. 


= pulse repetition frequency. 


Sh 
ot 
' 


n 
iz) 


= scanning frequency. 


> 
be 
" 


bandwidth for noise purposes. 


~ 


= input circuit bandwidth. 


Cc Pb 
SH 
H 


pear combined R F and I F bandwidth of continuous-wave-system receiver. 


ig = bandwidth multiplying factor = 1 for simple L C circuit. 
¥y = number of pulses received during detection time. 

yi = y/N 

Be ace mutual conductance of first receiver tube. 


= gain of transmitting antenna. 


= height of radar antenna. 


x 


= target height. 


~~ 


modified Bessel function of the first kind. 
= Boltzmann’s constant. 


= wave length of transmitter. 


Saher he We an wi a 
~ 
N 
~~ 
iT] 


= sweep length in miles. 


Pee teed 587 

n' =n/N 

” = number of pulse intervals per sweep. 

7 = number of separate velocity channels in continuous-wave-system receiver. 

N = number of pulses integrated, or, in cw system, the number of independent vari- 
ates (of length Webs.) integrated. 

We = number of pulses per scan. 

NF = overall noise figure of the receiver. 
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Y 


Py = probability that N pulses of noise will exceed a given level. 

P = probability that N pulses of signal plus noise will exceed the bias level. 
P' = probability that at least one group of N integrated pulses will exceed the bias 
level within the detection time. 

Ps = average power. 

Lig = transmitted power. 

Pin ~ minimum detectable power at receiver. 

p. = effective input noise power to receiver. 

r = resistance. 

R = radar range. 

Riax = Maximum radar range. 

R, = idealized radar range. 

Re = equivalent noise resistance of first receiver stage. 

R, = total shunt resistance of first receiver input circuit. 

AR = range interval for integration with a moving target. 

o = scattering cross-sectional area of target. 

9 ier pulse length. 

tee false alarm interval. 


a4, 5° detection time. 

ae maximum integration time for moving target. 

T = absolute temperature. 

7, = absolute temperature of space radiation :received by antenna. 
Tp = absolute temperature of room. 

v = velocity of the target. 


velocity of traveling gate. 


S 
) 
i] 


= visability factor of receiver. 


angular velocity of antenna. 


th RS as 
wv 


= received signal pulse energy in units of k Tp NF. 


-noise level in units of the rms value of noise — the bias level. 


= 
" 
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A STATISTICAL THEORY OF TARGET DETECTION 
BY PULSED RADAR 


PART I - INTRODUCTION 


THE USUAL RADAR RANGE EQUATION 


Most radar engineers are now well acquainted with the following equation used 
to determine the maximum range of a pulsed radar system: 


P GA oa $ 

R 1 t e 

max IP 16 7? (1) 
min « 


where | 
= peak transmitted power in watts. 


= minimum peak detectable signal in watts. 


o = scattering cross section of target in units consistent with range. 

G = gain of transmitting antenna. 

A effective area of antenna for receiving in units consistent with 
range (usually about 2/3 of the physical aperture, A, = GM /47), 

a dimensionless loss factor which accounts for atmospheric absorption, 

_ losses in antenna and transmission lines, etc. 


om 
" 


The number of pitfalls that may be encountered in the use of the above equation 
are almost without limit,’ and many of these difficulties have been recognized in the 
past )» 8) Three of the most troublesome are: 


1. The Scattering Cross Section 


In the case of moving targets, the wide variation of this quantity with as- 
pect, and hence with time, is a matter of vital concern. The variation of. 
cross section as a function of frequency may also be critical. 
2. The Minimum Detectable Signal, P,. . 
The statistical nature of the noise with which Pe must compete makes this 
an ill-defined quantity. 


' For references see page 143. 
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3. The Maximum Range, R 


max ‘ 


The statistical nature of P|;| in turn makes Loa a statistical quantity. 


There are also lesser troubles, such as the dependence of 8, the loss constant 
on the range, and the contribution of reflections from the ground, sea, or other ob- 
jects to the incident and received powers. (One must also remember that a target 
cannot ordinarily be detected at ranges (in miles) much greater than V2h_ +N2h,, 
where h_ is the height of the radar antenna and h, the height of the target in feet, 
except in the case of superrefraction, or "ducts." See pp. 55-58, Ref.(18). This is 
the familiar "line of sight" limitation due to the earth’s curvature. 


THE SCATTERING CROSS SECTION OF THE TARGET 


For a stationary radar observing a stationary target, the scattering cross sec- 
tion is a constant. Although it may not.be calculated for any but the most simple 
target shapes, it is not too difficult to measure. On the other hand, if either the 
radar or the target is in motion, the cross section becomes a function of time caus- 
ing the return signal strength to fluctuate. In general, the plot of cross section 
as a function of angle for a complex target such as an aircraft shows two interest- 
ing features. There is a nearly continuous rapid fluctuation having an angular pe- 
riod in the neighborhood of a degree or so (for \ in the microwave region), and a 
slow variation with a period in the order of 20° or more. Both of these variations 
may be as great as 30 db. The question at once arises: In lieu of using the com- 
plete polar diagram of cross section vs. angle, what kind of average figure can be 
used, and under what conditions? The answer to this question involves such things 
as angular rates of the aircraft with respect to the radar, correlation times, rep- 
etition rate of the radar, and number of pulses integrated. It is almost obvious 
that the only general way to treat this complex problem is to consider the cross sec- 
tion as a statistical variable. This approach seems mathematically feasible. However, 
in the present report the cross section will be considered to be a constant. An at- 
tempt to justify this assumption is the following: The rapidly fluctuating corre- 
lation angle at half-power points is perhaps 1.1°. The normal variation in attitude 
angle of an aircraft may be about 30° per second. (This variation may be caused by 
small rapid changes in pitch or roll due to normal turbulence of the air as well as 
by systematic changes in position.) Thus, the corresponding correlation time for o 
is around 1/300 second. If the observation time is essentially greater than this 
period, it may be assumed, as a first approximation, that the rapid fluctuations in 
the cross section "average out." 


The slow variations (period around 20°) may or may not average out. However, if 
the average over all likely attitudes is used for 0, or tobe more exact, if a weighted 
average is taken for o according to the probability for any attitude, then the prob- 
ability of detection may not be changed very much. Henceforth, in this report o 
will be assumed to be a constant, on the basis of the above statements. It may be 
mentioned in passing that o loses its meaning if the target is not uniformly illum- 
inated. Such can be the case, for example, if waves reaching the target via two or 
more paths combine to produce an interference pattern at the target. This effect 
exists in the detection of ships by surface radar. 
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THE MINIMUM DETECTABLE SIGNAL 


As is well known, ‘1»?»%) the minimum detectable signal power in aradar receiver 
is fundamentally limited by three main factors; i.e., Johnson noise in circuit ele- 
ments of the input circuits, shot effect and other noise in the first tube (and to 
some smal] extent succeeding tubes), and cosmic noise picked upby the antenna. There 
may also be man-made interference such as engine noise, radiations from other radars 
and radio transmitters,etc Clutter caused by sea return, rain, clouds, land masses, 
etc., may reduce the minimum detectable signal by a considerable amount. The effects 
of clutter and man-made interference are complete subjects in themselves, ‘'®) and 
will not be treated further in this paper. A study will be made here of radar range 
in the absence of such interference. It is not too optimistic to suppose that cir- 
cuits will eventually be designed which will largely eliminate man-made interference, 
and most types of clutter. 


The mean squared noise voltage across a resistor of resistance r is given by 


e? = 4kTr A f (2) 
whe re 
k =Boltzmann’s constant = 1.38x107?° joules/degree 
T = absolute temperature of the resistor 
Af = the frequency interval under consideration. 


Though the noise at the input circuit of a receiver is usually several times this 
value, it provides a convenient scale for measuring the input noise. The effective 
input noise power is defined to be 


p = kT, Af NF (3) 


where NF is the so-called noise figure of the receiver, and Tp is the absolute room 
temperature.* Ifa signal power of the same value as p were incident on the antenna 
and the receiver were noiseless, then the output would be the same as in the case 
when noise only was present. 


: ; ; 5 : 3 1 
At this point, one important result concerning the noise figure due to Herold‘? 
is pertinent: 


NE eee ot ae f(R,) (4) 


Le} 
ti] 


absolute temperature of space radiation received by the antenna. 


<2 
" 


room temperature. 


* Complete discussions and derivations will be found in the Mathematical Appendix (a separate 
report), ‘* 
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> 
— 
" 


bandwidth of input circuit. 


‘ 


C = total shunt capacity of input circuit. 

Reg = equivalent noise resistance at input (due mainly to shot noise in 
first tube) © 2.5/g_ for triodes. 

F =a factor depending on the exact type of input circuit coupling (= 1 
for simple tuned circuit). 

R = input shunt resistance including effect of finite input resistance 


of tube. 


f(R,) =a function of R,, R C and Af’. 


This formula assumes a more or less conventional type of input tubes, such as the 
VHF triodes and pentodes. However, it seems reasonable to believe that the gen- 
eral conclusions which are reached from Eq.(4) will apply to velocity-modulated input 
tubes as well. 


The main points to be noted about Eq.(4) are these: 


it f(R,) approaches zero as R, approaches infinity. R, may be increased by bet- 
ter tube design. 


2.C/g, should be made as small as possible inatube used as the first amplifier. 


3. Long pulses tend to allow smaller bandwidths for the input circuit, and hence 
lower noise figures. 


4, If R, CAf' is made small enough, and R, large enough, the noise figure will 
approach T/T. r 


Point 4 is of the greatest importance. It sets alimit on the noise figure when there 
are no sources of noise in the receiver itself. Though such a receiver will never be 
built in practice, it may be possible some day to approach closely this ideal state. 
Then the input noise will be almost entirely dependent on the temperature of space; * 
or, in other words, on the noise received by the antenna from without the radar set. 
That this state of affairs is not yet at hand is evidenced by the fact that at present 
the noise figure for microwave receivers is around 10, and for longer-wave receivers 
perhaps as low as 3 or 4. 


The concept, often stated, that the ideal noise figure of a receiver is 1.0 is 
erroneous.** This would be true only if the temperature of space were the same as 


room temperature. Actually the temperature of space decreases rapidly with decreas- 
ing wave length. ‘) 


* Though the noise figure can be decreased by increasing T,, this would increase the actual 
input noise, as is apparent from Eq. (3). 


** The noise figure of a receiver may be defined in such a way that the antenna must be re- 
placed by a resistor at room temperature equal to the radiation resistance of the antenna. 
In this case the ideal noise figure of the receiver would be 1.0. 
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The average space temperature* is around room temperature at 180 megacycles and 
drops to around 30° absolute at 450 megacycles. ‘1?) No good measurements are avail- 
able in the microwave region, but there is reason to believe that values of 10° or 
lower may be found. If this proves to be true, then it is conceivable that the noise 
figure of future microwave receivers may be improved by a factor of 100, which would 
mean that the range of radar sets could be more than tripled as a consequence of this 
one factor. It is certainly a field where research should be pushed to the utmost. 


It has often been the practice to calculate the maximum range of a radar set 
from (1) by assuming that pee = kTp Af NF, or that the minimum detectable signal power 
is just equal to the average noise power. 


This gives y 
PGA oS 
max 16r? kT, A fNF 
Now the energy per pulse is represented "by 
E. = Lead (6) 
where un is the pulse length. Making this substitution in (5) gives 
% 
ae Pime 
Snax | JeontaT AF | G7 ap® 
WF | Gap) 


It is usually said that if 7 Af is made equal to 1, the amplitude of the pulse after 
passing through the amplifier will not differ much from the amplitude which would re- 
sult if the pulse were infinitely long. Without further ado, 7 Af is put equal to 
1, and the resultant equation 


: y% 
E GA 08 
R = eee ee (8) 
max 167° kT, NF 


emerges as the radar range equation. Now the unfortunate fact (in some respects) is 
that the range of a radar set calculated by means of this formula often turns out to 
be rather close to the experimental range. Naturally, under these circumstances great 
effort has not been expended in investigating the validity of radar-range equations. 


* There is a variation of the space temperature with direction(!2), When the antenna points 
near the horizon, the temperature may be higher than when it is pointed at the zenith. In 
particular, if any appreciable part of the radiation strikes the ground, the thermal radi- 
ation received from those directions will have a temperature nearly equal to the actual 
temperature of the surroundings. 
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The reasons for the agreement of equation (8) with experiment are many. First 
of all, the cross section has been, in most cases, determined by observing the max- 
imum range of a particular target and solving equation (8) foro. This one fact 
alone accounts in no small way-for the agreement. Secondly, the fourth power law 
makes R rather insensitive to changes in the various parameters concerned in equa- 
tion (8). A much fairer test is to compare respective values of | eg rather than 


R 


max 


Equation (8) is in no sense perfect with regard to its agreement with experiment. 
Errors of as much as + 30% are common, and factors .of 2 can often be found. However, 
considering all the unknown factors present in an experimental determination of max- 
imum range with an operational radar set, this agreement is considered to be quite 


good. 


In any field of science, theoretical equations are deduced to explain observed 
data. However, one is very cautious in using these equations to predict results for 
other experiments where the values of many of the variables differ greatly from those 
used in the particular experiments already performed. Most of the radar sets built 
to date have operated within essentially narrow limits as far as some of the para- 
meters are concerned. Particular examples are pulse repetition frequency, and, most 
important, the number of pulses integrated. This latter quantity is not even men- 


tioned in equation (8); but, as will be seen in the next section, it is of vital im- 
portance. 


The task is now two-fold: 
1, To make a satisfactory statistical definition of the range of a radar system. 


2. To determine the dependence of this quantity on the parameters of a (pulsed) 
radar system. 
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PART II 


THE STATISTICAL PROBLEM OF THE MINIMUM DETECTABLE SIGNAL 
AND THE MAXIMUM RANGE 


GENERAL BACKGROUND 


It has been realized by many workers in the field* that the range of a radar set 
is a statistical variable and must be stated in terms of probabilities rather than in 
the exact terms of an equation such as (8). However, the evolution of a practical 
working theory does not seem to have been accomplished so far. The following work is 
a first step in that direction. 


Before beginning the explanation of equations and derivations, it will be well 
to glance at some of the new ideas which will be included. 


The random noise, which limits the range, can at intervals assume large values 
due to its statistical nature. This means that there will inevitably be times when a 
random fluctuation of the noise will be mistaken for a signal. The average interval 
at which such undesirable events take place will be called the false alarm time, and 
it will be found that the probability of detecting a target will be a function of 
this time. Let the reader at once be cautioned ‘against thinking, “If it were a noise 
flash, I can easily tell by looking a little later. If it were a signal, it will 
still be there; if it were noise, it will be gone." 


The second new parameter which will be introduced is the detection time. It is 
apparent that if an observer can spend sufficient time in deciding whether or not a 
target is present on an oscilloscope screen, the probability of a correct decision 
being reached will be increased. It is also obvious that in any practical situation 
in which radar is used one cannot take unlimited time to decide whether or not a tar- 
get is present. To put things on a quantitative basis, the time in which a decision 
shall be rendered must be specified. In this event, there will not always be time 
for the "second look" just mentioned; but should time permit, then the probability of 
detecting a target will be increased at the expense of a longer detection time. Even 
so, there will still be a certain lesser probability that the noise flashes will oc- 
cur on both occasions. Further, it will be found that the velocity of a moving tar- 
get has an appreciable effect on the detection probability, due to the fact that the 
signal from such a target does not "remain stationary" (see page 80). 


* For an excellent qualitative statement of the problem, see Radiation Laboratory Series No.l. 
pp.35-47, Ref.(18). 
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PRELIMINARY STEPS ‘ 


It is desirable to present data in the most compact form, and the first step in 
this direction is the elimination of the necessity for the appearance of such para- 
meters as E_, G, A en On $, and NF in the final results. To this end, a parameter R, 
is defined ehich is given by a alight modification of Eq. (7), as follows: 


E GA o8V | ™ EGA. o8V 7 
mm | eee Ree Rhee ae (9) 
16r°E, : 1677*kT, NF : 


Here, the factor 1/7, Af has been replaced by V, the so-called visibility factor.* 
This factor will always be less than 1 but usually not less than 0.8, except when the 
Doppler effect is very large. WR, will be called the "idealized range" for lack of a 
better term. 


Now let the received energy per pulse at any range R be E,. Then it is clear 
from the equations (9) that 


% 
kT,NF 
ave ee se (20) 
Ry Ep 
and defining 
E 
pe (11) 
k R 
gives from (10) 
Be ee 
R, x 12) 


* The derivation of exact formulas and numerous curves of visibility factors as a function of 
pulse width, bandwidth, type of amplifier, and off-resonance of carrier frequency will be 
found in the Mathematical Appendix (a separate report) (23), The visibility factor is actu- 
ally given by 


Raf 


where E., is the maximum voltage to which the pulse rises at the receiver output, and E, 


2 
is the steady state voltage at the same point. The quantity ( E should be contained 


ss 
in (7) and (8) but is usually omitted because it is so near to unity when 7Af=1. In the 


case where the bandpass characteristic of the, pape is the conjugate transform of the 
pulse, the visibility factor is exactly unity 
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where x is now the signal pulse energy in units of the average receiver noise pulse 
energy. As an example, suppose x=4, which means that the signal power equals four 
times the average noise power. Suppose the probability is calculated to be 0.5 that 
in this case the signal will be detectable. There is then a point P=0.5 at R=0.7R. 
When a series of such points are calculated for various values of x, a curve for P e 
a function of R/R, may be drawn, assuming fixed false alarm time, etc. 


INTEGRATION OF PULSES 


Before proceeding further, the meaning of pulse integration must be defined in 
detail. In its simplest form, it merely consists of adding N successive signal pulses 
together and attempting to detect the sum rather than an individual pulse. Now, what- 
ever the integrating device may be, it will not know in advance whether there is a 
signal or not, and hence in the absence of a signal it will add up N successive noise 
pulses. Therefore, the comparison is between N signal plus noise pulses and N noise 
pulses as contrasted to a single signal pulse to a single noise pulse. One might be 
tempted to say that the signal to noise ratio would be unchanged, and that integra- 
tion, or addition, of pulses therefore offered no advantage. This argument neglects 
the fact that the noise voltage fluctuates about its average value The mean or aver- 
age value of the noise voltage is not of too much concern, for it can always be "bi- 
ased out.'"* If we add N signal pulses of voltage V, the total signal voltage is NV. 
If we add N noise pulses of average voltage Vy, the average sum will be NV,. How- 
ever, the average sum can be balanced out. The question is, whether or not the fluc- 
tuation in the sum voltage is now N times the fluctuation voltage of single pulse. 
Tf the answer were yes, then integration would be futile. However, due to the random 
nature of the fluctuation of any single pulse, the fluctuation voltage of the sum is 
only about VN times the fluctuation voltage of a single pulse. It is the signal to 
noise-fluctuation** ratio, not the signal to average noise ratio that is of paramount 
importance. The greater the number of pulses integrated, the greater is the signal 
to fluctuation ratio, and the greater is the probability of detecting the signal, but 
at the expense of longer detection times. 


DEFINITION OF DETECTION AND THE BIAS LEVEL 


Before the false alarm time can be calculated, a definition of "detection of a 
signal" must be given. Detection of a signal is said to occur whenever the output 
of the receiver exceeds a certain predetermined value hereafter called the bias level. 
In the absence of any signal, this bias level will on occasion be exceeded by the 
noise alone. The higher the bias level is set, the more infrequently this happens. 
The first problem is to calculate the required bias level, given the false alarmtime. 
Knowing ‘this bias level, the rest of the problem is to calculate the probability that 
any given value of signal (plus noise) will exceed this level. 


* Practically, the bias level should not be too large, or the fluctuations in the bias will 
become of concern. See page 78\where a method of reducing the necessary bias level by a 
considerable factor is discussed. 


** The mathematical term for the fluctuation is the "standard deviation," usually denoted by o 
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This is well and good, one says, but is thiS the best means of detection? What 
about the operator watching a cathode ray tube - what are his criteria for calling 
"signal"? Of course, it is impossible to say exactly, as is evidenced by the wide 
variation among radar operators. One can see, though, how an operator is affected 
by the false alarm time. If he is told that he will be subject to severe penalties 
if he calls a false alarm (calls a signal when it subsequently turns out that there 
was none), then he will be very cautious. If a doubtful pip appears on the screen, 
he will use discretion and say nothing. This means that under these conditions the 
false alarm time is increased, and at the same time the probability of detecting a 
target at a given range is decreased. 


The operator may use the shape of a signal pulse contrasted to that of the noise 
as a criteria for detection as well as amplitude differences. This is thought to be 
a second order effect. The operator, on the other hand, is limited to some extent 
by the minimum brightness ratio which the eye can detect. 


It seems that the method of electronic detection proposed above will be practi- 
cally as good as any other possible method, electronic, human, or otherwise, if 
identical false alarm times and detection times are assumed. This statement is cer- 
tainly not to be considered obvious. It should be possible to make some experiments 
to verify this theory. 


METHODS OF PULSE INTEGRATION 


As stated before, to integrate pulses it is merely necessary to add them to- 
gether. There are many different practical ways in which this is done. One of the 
simplest is the use of a cathode ray tube screen‘®), Due to the screen persistence 
time, a certain number of pulses will be effectively integrated. In this case it 
will not be a simple addition, being more in the nature of a weighted average. The 
effect of weighting is always bad. In other words, the effect of equal samples in 
the integrated result should be as nearly the same as possible. PPI _ type of pre- 
sentations which use intensity modulated displays usually have much longer integra- 
tion times than an A scope. 


One must not overlook the human operator,* who goes along with the cathode ray 
tube, as a vital part of the detection mechanism. The combination of the eye and the 
brain makes avery good integrator. In fact, the maximum integration time for a skil- 
led operator may easily be several seconds. The best electronic integrators for pulsed 
radar built to date wil] not better this figure to any great extent. Henceforth, a 
mode] electronic integrator which linearly adds N pulses will be assumed. 


* There are a large number of factors involving observers and oscilloscopes which are quite 
complicated and are more or less outside the intended scope of this report. Lawson and his 
group have done a great deal of work on this subject, the results of which will appear in 
Chap.VIII of Ref.(19). Most of these experimental results are also available in Ref.(24). 
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Now, pulses can be integrated in the RF stages, in the I F stages, or in the 
video stages (1), (14), (16), (18) Futhermore, there canbe one or more linear or square 
law detectors present*, and the integration can be done in one or more steps and in at 
least two different ways. Manyof these possibilities are reserved for detailed treat- 
ment in a separate mathematical report (#9). 


Fortunately, the results for the various cases show little difference, with one 
marked exception. RF and IF integration are better than video integration for 
small signals (compared to the noise). However, there is no practical way known at 
present to take full advantage of RF or I F integration with moving targets be- 
cause of the requirement that the successive received pulses must be completely co- 
herent ‘*9)» (1*) | Coherent integration would be possible in the case in which both 
the radar and the target were stationary, but this case is not of much practical value. 
The difference between various types of both detectors and video integrating circuits** 
is small, as far as results of this kind of study are concerned. There are, of course 
Many reasons why a choice is made in practice, such as sensitivity to small changes 
in amplifier gains, vulnerability to countermeasures, etc. 


It is worth describing one scheme for integrating in which a pulse known to be 
only noise is subtracted from each possible signal plus noise pulse. WN of these com- 
posite pulses are then integrated. With no signal, the average value of any number 
of such composite pulses is nearly zero, so that the required bias level is consid- 
erably reduced. Such a method is much less sensitive to a small change in bias level, 
and would usually be preferred in practice. This case is much more difficult to 
calculate than the straight addition case; and since sample calculations show the 
results to be nearly identical, the latter method has been used to obtain the curves 
of Figures 1] thru 50. 


Figures 51 and 52 show the difference in sensitivity to bias level for this 
method. Figures 53 and 54 show the comparison of straight integration to the case 
in which a noise pulse is subtracted from each signal-plus-noise pulse. 


Practical types of electronic pulse integrators often take the form of very nar- 
row band audio filters having-their center frequency at the pulse repetition fre- 
quency ‘?1) or some harmonic thereof. The action of such a filter can be understood 
roughly by consideration of the frequency spectrum of a finite group of N pulses. The 


* It is assumed throughout this report that the video bandwidth is large compared withthe I F 
bandwidth. Actually, the results will be affected only if the video bandwidth is small com- 
pared with the I F bandwidth (24), a condition not often found in practice. 


** One might ask if there would be any advantage in having an integrator which adds the sum of 
the squares of the N pulses or perhaps the sum of some other function of the amplitude. Ac- 
tually, it can be easily shown that this just corresponds to changing the shape of the de- 
tector curve, and what is being asked is, "Is any shape of detector curve much superior to 
the linear or square law form?" Apparently the answer is no. There is a "best" detector 
curve. for every different signal strength, x, given by log Ig(v v2x) where I, is a modi fied 
Bessel function. No results have been obtained for this detector function, but it is thought 
that the maximum difference in range between this and the square law or linear detector will 
not exceed five percent. 
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envelope of such a spectrum is simply the familiar sin x/x curve of a single pulse, 
while the actual curve has appreciable values only in the neighborhood of the harmon- 
ics of the repetition frequency (including dc)*. The greater N is, the more closely 
the spectrum clusters around these harmonics. Thus, the filter may be made narrower, 
excluding more and more noise, but retaining most of the signal energy. With such a 
narrow band filter-type of integrator it is very simple to subtract a noise pulse 
from each signal-plus-noise pulse by gating the receiver at a frequency double the 
center frequency of the filter. To prevent the possibility of a signal on every other 
gate, the sweep length would ordinarily be held at less than one-half of the pulse 
repetition period. The simple electronic type of integrator has the disadvantage of a 
fixed integration time. If the number of pulses returned from a target is greater or 
less than the number of pulses for which the integrator is set, the operation suffers. 
With the human operator, the story is different. He can adjust his integration time 
rapidly to fit changing situations. This procedure could be approximated electroni- 
cally by the use of two or more successive integrators in series,** or by the use of 
so-called "weighting circuits." Such a complicated procedure does not come within 
the scope of this report. 


METHOD OF OBTAINING THE BIAS LEVEL 


By means covered in detail‘?*), in a separate mathematical report***, the probabi- 
lity that the sum of N pulses of noise voltage alone will be greater than an arbi- 
trary level y is obtained. This relation may be symbolically represented by 


Pe ry) (13) 


where y is measured in units of the rms value of the noise. The number of groups of 
noise pulses which are observed in a fixed false alarm time, Ty,» 1s then found. 


Wher speaking of noise pulses, it is convenient to assume mentally a range gate 
equal to the pulse length at a fixed range. If the range sweep is continuous, such 
as with an A scope, the effective number of independent noise pulses observed in one 


* There is a close resemblance between such a spectrum and the diffraction pattern of an N 
slit grating (see any standard text book on physical optics). 


** The advantage of a multistage integrator is that if a signal which is large enough so that 
the number of pulses which need to be integrated in order to produce a detectable signal 
occur in a time appreciably less than the total integration time, one of the sub-stages 
will detect the signal much sooner than will the final stage. 


*** Tt turns out that the functions which describe the probability that the noise alone, or a 
given strength signal plus noise will have any arbitrary amplitude, are quite complicated 
and hence only some of the results and general procedures are given in this report. Fur- 
thermore, it should be mentioned in passing that the use of the central limit theorem, or 
the so-called "normal approximation," is not valid until the number of pulses integrated is 
of the order of 1000. This is because the values of the distribution functions far out on 
the tails play a major role in the calculations. Several investigators in the past have 
made the mistake of assuming that the normal. approximation was satisfactory if N were only 
of the order of 10. 
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repetition period is given by the length of the sweep divided by the pulse length* 
hereafter called 7. It is apparent that 7 = 2L/cTp = 10.8L/r, where L is the sweep 
length in miles, c is the velocity of light, and 7. is the pulse length in micro- 
seconds. In the special case inwhich the sweep occupies the total time between pulses, 
he V/T fh which is merely the reciprocal of the duty cycle. The time for N pulses 


e ap 
é fers tofesr : 5 ; 
to occur is N/f Therefore N/f.~  N  8roups are observed in the time T pq? assuming 
r 


only one gate per sweep. Since the effective number of gates per sweep is 7, the to- 
tal number of independent chances for obtaining a false alarm in Fee is** 


Tak 
Ma ne (14) 


The false alarm time is defined as the time in which the probability is % that the 
noise will not exceed the bias level.*** From (13) and (14), 


Cigaipe eee (15) 
n 2 


from which y, the bias level, is obtained. 


PROBABILITY OF DETECTING A SIGNAL 


Having established the value of the bias level, the probability that a signal 
will exceed this level in a given time, namely the detection time Tq: must be calcu- 
lated. The signal is assumed to consist of N integrated pulses. The time of such a 
pulse group is N/f The number of such groups which occur in iP is given by 


(16) 


As a corollary to the previous definition of detection, it is now assumed that the 
signal is detected if any one of the y’ groups of pulses exceeds the bias level. One 
will ask, at this stage, "Why not count exactly how many times the signal exceeds the 


* If the range gate is much wider than the pulse length, the operation of the integrator will 
suffer more or less, depending on the exact type of integrator used. This corresponds 
somewhat to the case of an oscilloscope where the spot does not move by at least its 
diameter within a pulse length. 


** This derivation assumes that the antenna is not scanning. With a scanning antenna, inte- 
grating channels must be deposed in angular position as well as in time. In order for 
(14) to hold, the number of pulses per channel per scan must be equal to or greater than 
N, the number of pulses which each channel integrates. 


*** This very nearly, though not exactly, corresponds to the earlier definition given on p. 7]. 
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bias level?" This would in effect correspond to a two-stage integrator. Such a de- 
vice is not considered here, though it is easy to make an extension of the present 
theory to cover this case. 


At any range R, the normalized signal strength x is obtained from Eq.(12). The 
probability that the signal plus noise will exceed any value y for a single group of 
N integrated pulses is known (#9), and may be represented symbolically as 


P = f(y,x). (17) 


The probability that at least one of the y’ groups will exceed the bias level y is 
then 


P' = 1-(1-P)”. (18) 


Notice that yy’ must be an integer for the analysis to be strictly correct. It will 
be satisfactory, however, if one always requires y’ >1*. 


EFFECT OF ANTENNA SCANNING 


If the antenna is scanning, some modifications of Eq.(16) for y’, the number of 
groups of pulses integrated, will be necessary‘?*). If, with a P PI type of pre- 
sentation, the antenna moves at an angular velocity w, and the beam width is B, then 
the number of pulses per target per scan will be 


Bf 
renee D (19) 
sc (G3) 
and (16) is replaced by 
Tou 
yee Pe 5 eee SE 
9 MS Nie N (20) 


where f . is the number of scans per second. With a simple type of electronic inte- 
grator, N. must be equal to or greater than N for Eq.(20) to be valid, assuming that 
the integrator does not hold over from scan to scan. If the integrator does hold over 
from scan to scan, as an operator partially does, then it is only necessary to have 
y'21 as before. In any case (20) only holds if Taf cael 


VG Tea, then y’ = ie /N, which must be equal to or greater than 1.* 


* It is always best for y' to equal 1. In this case the integrator effectively integrates 
pulses during the whole of the detection time. ‘y>l is the case in which the detection time 
is longer than the integration time. Here the probability for detection is greater than if 
the detection time were reduced to the integration time, but less than it would be if the 
integration time were increased to equal the detection time. The case for y'<l is that one 
in which the number of signal pulses occuring are fewer than the number for which the inte- 
grator is set. In this case the probability of detection is reduced from the value it would 
have if the integrator were set for exactly the number of signal pulses which occur. Tocal- 
culate this latter case would require using N to calculate the bias level as in (15), but 
the use of some lesser value N in obtaining (17). This will be done, but results have not 
been obtained as yet. 
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PRESENTATION OF THE RESULTS 


The results are presented in the form of a set of curves. This is necessary he- 
cause of the complicated form of the analytical solutions. The parameters involved 
in the curves are: 


" 


P 
R/R, 


the probability of detecting a target at range R, 
the ratio of the range to the idealized range. 


i 


: - Grea Bast 


tae the false alarm time 
f, = the pulse repetition rate 
7 = the number of pulse intervals per sweep. 
= ° s * 
x Td dee Ne 


ey the detection time 


f,- 7 the scan frequency 


= 
" 


the number of pulses per scan 


= 
" 


the number of pulses integrated. 


A summary of the range of the variables for the curves presented will be found 
on page 84. 


AN EXAMPLE WITH A QUASI-STATIONARY TARGET 


A simple example is now solved assuming a stationary target. The radar set will 
also be assumed to be stationary. The following data are taken as given: 


w = angular rate of antenna = 30°/sec, ine = 1/12 
370° 


WW 


B = beam width of antenna 
f, = pulse repetition rate = $00 per second 
SS pulse length = 1 microsecond 
R 


> = idealized range for given target and average aspect = 40 miles 


T fa 
Td 


= required false alarm time = 5 minutes 


= required detection time = 25 seconds. 


elt ag fale if the sweep occupies the total time between pulses. 


** See also (16), and the conditions on (20). The notation used on Figs.1-54 is ¥ = LMA: 
which represents the special case in which there is no scanning. In general this should 


be replaced by Y= TBs ON eS 


7 


Type of detector - electronic integrator, ‘V = 50; sweep length = 20 to 80 miles. 


Step 1. Calculate N.. from (19) 


Bf, — 3x500 


Nia aay oe Baan 


Step 2. Calculate y from (20) 


1 — 
25x;5%50 = 104 


Step.3. Calculate 7 from 7 = 10. 8L/r, where L is the sweep length in miles and T is 
the pulse length in microseconds. 


10. 8x(80-20) 
eR BO a 


Step 4. Calculate n from (14) 


n= ee ec = (5x60)x500x648 = 0.98108 


Step 5. Refer to Fig. 23; n = 10® and y = 100. Mentally interpolate a curve for N = 
50 between N = 30 and N= 100. This curve gives probability of detection at 
any R/R.. R, is given as 40 miles. For instance, P = 0.50 at R/R, = 1.07 or 
at R = 43 miles. 


MOVING TARGETS AND/OR RADAR 


If there is an appreciable change of range with time between the radar and the 
target, a limit will ordinarily be set on the number of pulses which can be integrated. 
This is because the returned pulses will just fail to overlap when the target has 
moved through a distance d = 7_C/2 where c is the velocity of light. The effective 
distance over which the pulses can be assumed to contribute their full amplitude is 
about % this value. If the rate of change of range is v, the time available for in- 
tegration is 


7 (21) 


The maximum number of pulses which can be integrated in this time is 


Neo i Tey, 9 Av 3 (22) 
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This quantity Nax 1S the maximum number of pulses that can be integrated, provided 
that it is not greater than N,,: the number of pulses per scan. In the case where 
Nex e6? then Ni is the maximum number of pulses which can be integrated. 


In the case of approaching targets, one may be concerned with the probability 
that a target will be detected by the time it has reached a certain range. Assuming 
the target to have started its approach at range R,, the probability that it will have 
been detected at least once by the time it reaches range R is 


R 
p = 1-]] Pn) (23) 
R, 


where R progresses from R. to R in units of AR. The length of the AR intervals 


and the number of pulses integrated per interval are determined from the considera- 
tions given above. 


An example follows in which (23) can be reduced to a particularly simple form: 
Assume a continuously directed beam (no scan) and the target moving toward the radar 
with a constant range rate v. The finite product in (23) may be approximated by 


- 
: 1 

log, II [a-P¢R) | ap: S c-pyar (24) 

7 R, R, | 


and using AR = d =—7~ equation (23) becomes 


4 R 
7 s (1-P) dR 
P = l-e 1 (25) 


The integrations necessary in the solution of this type of problem must be performed 
numerically, using the graphical data of figures 1 to 50. 


In problems where the antenna is scanning, equation (23) may be approximated in 
different ways depending on the exact values of the parameters involved. ‘These are 
rather simple to work out in any specific case. 


* A system could presumably be built incorporating one or more velocity gates. Such a velocity 
gate would travel with a given preset velocity. In this case, the relative velocity of the 
target to that of the gate, v-vg , can be used in Eq. (22) in place of the target velocity, 
v. The greater the number of velocity gates used, the greater will the probability be that 
the difference between the target velocity and some one of the gates will be very small. 
Therefore, in this gate the allowable value of Nmax will be large, and the probability of 
detection in this gate will be increased. 


In any multi-channel receiver, such as this, the number of pulse intervals per sweep, 7, must 
be multiplied by the number of channels in calculating n. 
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EXACT EFFECT OF THE NUMBER OF PULSES INTEGRATED ON THE RANGE 


One might expect that for.a given n and a given probability of detection, the 
range to the fourth power would vary as N, as was stated on page 73. This would be 
true with coherent integration, but with video integration the variation is between 
N% and N (assuming a threshold signal). This effect is due to the so-called "modula- 
tion suppression" of the weak signal by the stronger noise in the process of detection. 


Fig. 55 shows the exact variation of the exponent of N, here called @, as a 
function of N, and of n, for P fixed at 0.50. The effect of n is seen to be quite 


small. 


Fig. 56 shows the variation of the exponent of N for an incremental change of N 
as a function of N and n. P is again fixed at 0.50. In both cases, 6 approaches 0.5 
as N approaches infinity; though much more slowly, in the first case. 


APPLICATION OF RESULTS TO CONTINUOUS-WAVE SYSTEM 


Though this report is concerned primarily with pulsed systems, the results are 
directly applicable to continuous-wave systems. To accomplish this, the following new 
notation is introduced: 


= the average cw transmitter power. 
Af. = the combined RF and I F bandwidth of the cw receiver. 


: the number of separate velocity channels incorporated in the receiver.* 


3 
t 


The change-over is then made by means of these substitutions: 


Replace E, by Po Aft » in R 


Paci 
Put n = Tg? OS og 


N is now to be taken as the number of variates (of duration 1]/ Af.) which are inte- 
grated after detection. ** 


* In both the pulse and cw analysis it has been assumed that the range or velocity gates or 
channels are fixed in position. In the case where such gates’ sweep as a function of time 
in order to conserve apparatus (or for any other reason), the analysis is not strictly 
valid. A good rule-of-thumb is that the gate should move through the amplifier pass band 
in a time equal to the reciprocal of the amplifier pass band. In this case the effective 
visibility factor is about 0.8. Curves of the visibility factor for other sweep speeds are 
given in the Mathematica] Appendix (29) (a separate report). 


** Integration of M variates before detection merely corresponds to narrowing the RF (or I F) 


bandwidth by a factor of. 
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N must be less than y for the theory to hold. An optimum cw system is one in 
which y = 1 (7, = VAISS Ba ie and N= 1]. This gives the greatest range for a given 
energy expended during the detection time, 7, This corresponds exactly to the case 
N = 1 and y = 1 in a pulsed system. If the number of range channels 7 in the pulsed 
system is equal to the number of velocity channels 7’ in the cw system, then the two 
systems, with N = 1 and y = 1, will have identical ranges for the same average power. 


In either case, if N>1, a larger amount of average power is required, every- 
thing else remaining equal. In the pulsed case, reducing N necessitates higher peak 
powers, which may be impracticable; or it necessitates longer pulse lengths, which 
reduces possible range-resolution and at the same time aggravates the effect of a 
fixed Doppler shift due to the narrowing of the receiver pass band. In the cw case, 
reducing N necessitates a target with reasonably constant velocity so that the signal 
will not wander in and out of the pass band of the receiver, and also a sufficiently 
slow scan so that each target "pulse" is at least as long as the reciprocal of the 
receiver pass band. 
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RANGE OF VARIABLES FOR FIGURES 1 THRU 50 
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ERRATA 


In Figs. 13 thru 24, all of the ordinates appear as 
percentages but are labeled as probabilities. Therefore, 
in order to make the two scales conform, the decimal 
place should be moved two units to the left on all the 
ordinates. 
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SUMMARY 


In a previous report ‘?®) a statistical theory of radar 
detection was presented in outline form. The mathematical 
details were omitted, in order that the main ideas and results 
might be made available as soon as possible. 


This report contains the mathematics that led to the 
results presented in Ref.28. 


In addition, several subjects are briefly discussed 
that were not covered in Ref.28. These are collapsing loss, 
antenna beam shape loss, the effect of signal injection, limiting 
loss, ang moving target indication. 


ED, 


For references see page 264- 
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a 


n 
n 
N 


SYMBOLS 


amplitude of sine wave relative, to R.M.S. noise level 
one of the independent variables in the function Q(a, A) 
i*™ central standard moment 

one of the independent variables in the function Q(a, 8) 

half-power antenna beamwidth 

coefficient in the Gram-Charlier series 

characteristic function 

delta function 

base of natural logarithms 

frequency 

confluent hypergeometric function 

Campbell and Foster notation for characteristic function 

Campbell and Foster notation for anticharacteristic function 

the gamma function 

probability that the sum of N noise variates will exceed the bias level 
i*® Hermite polynomial ) 

index, subscript, or V-1 

incomplete gamma function as defined by Pearson (®) 

modified Bessel function of the first kind 

Bessel function of the first kind 


, th 


t cumulant 


standard i*® cumulant, or sometimes a modified Bessel function of the 
second kind 


integration loss 

collapsing loss 

generalized Laguerre polynomial 

number of excess noise variates integrated with N signal plus noise variates 
false alarm number 

n/N 


number of variates integrated 


*This symbol has a different meaning in RA-1506l. 


**This symbol is used in more than one sense in various places, but other meanings should 
be obvious. 
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p* = tw= wv-] ; 
P = sine wave amplitude 


pees probability 


Ps = probability that noise will exceed the bias level at least once within false 
alarm time , 
Py = probability that the sum of N ‘variates of signal plus noise will exceed the 
bias level 
pt = it® derivative of the error function 
Yo = R.M.S. noise level 


Qa,6)= modified Lommel’s function 


R** = envelope amplitude or radar range in R/R 


) 
Jay = idealized radar range ; 
fe) = collapsing ratio, ratio of total number of variates integrated to those 
containing signal 
s = cathode ray writing speed 
o* = standard deviation 
7 Cis (y-y)/o, semi-independent variable in Gram-Charlier series 
T,* = incomplete Toronto function 
M; = it® moment about the mean 
U. = Lommel’s function 
v* = normalized envelope amplitude 
v = tt moment 


w(f) = power spectrum 


w = 2nf 

C4 = power signal-to-noise ratio 

y* = normalized detector output 

y** = integrator output for the sum of N variates 


Y, = bias level 


*This symbol has a different meaning in RA-15061. 


-**This symbol is used in more than one sense in various places, but other meanings should 
be. obvious. 
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A STATISTICAL THEORY OF TARGET DETECTION 
BY PULSED RADAR: MATHEMATICAL APPENDIX 


BASIC FORMULAE RELATING TO THERMAL NOISE 


Both the thermal noise voltage across a resistor and the noise voltage due to 
the shot effect in a vacuum tube approach a normal distribution when the number 
of electrons involved per second in the processes tends toward infinity. In prac- 
tice, it may usually be assumed that the total noise voltage between any two points 
due to any combination of thermal, shot, and cosmic noise sources can be represented 
by the distribution function 


V 

Ly oa 
: dV l 
dP aris (1) 


where Yo is the mean square value of the noise voltage(*®), This distribution is 
valid provided all elements involved in the composition of the total noise voltage 
have been linear. 


If such noise is now passed through a linear filter of center frequency f, 
having a pass band which is narrow compared to f, the output will have an enve- 
lope, which has a probability density function 


R : 
cepa is (2) 
dP Yo e dR 


For references see page 264 
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where R is the amplitude of the envelope and ¥, is the mean square noise voltage, 


given by 
Y= i) w (f)df (3) 
F 


w(f) is the so-called power spectrum of the filter and is simply the square of the 
absolute value of the amplitude transfer function of the filter. 


If the input to a filter consists of a sine wave of frequency f,, as well as 
noise, then the probability density function of the output envelope is* 


R*+p2 
ahs e 
I, 


R 
dP = joe 7)eR R>0 (4) 


dP = 0 R<0 


where P is the amplitude that the sine wave would have at the output of the filter 
in the absence of noise, and J, is a modified Bessel function of the first kind 
(see footnote, page 167). 


The envelope of the output has a correlation time which 1s approximately equal 
to the reciprocal of the bandwidth of the filter. In simple language, it is im- 
probable that the envelope will change by an appreciable percentage in times much 
less than the correlation time, but it is quite probable that it will change by a 
goodly percentage in times large compared with the correlation time. It is probably 
a good approximation to assume that values of the envelope 1/Af seconds apart are 
independent, where Af is the bandwidth of the filter. By assuming such a discrete 
process it is possible to materially simplify calculations which would be very 
tedious if exact integration processes were used, while at the same time sufficient 
accuracy is obtained for most practical purposes. 


A further justification for this assumption in the pulsed case shows in the 
results. Changing the factor 1/Af to k/Af for the correlation time has only the 
effect of changing the false alarm time by the factor k. The probability of detec- 
tion turns out to be a very insensitive function of the false alarm time, so that 
if k is any factor of the order of magnitude of unity, the results are affected to 
a negligible extent. 


* It is of some interest to note that the same form of distribution function occurs in other 
problems. For instance, if yj represents the mean square velocity of a gas due to ordinary 

_ turbulence, and P represents the translational velocity of the whole mass of gas relative 
to some fixed reference, then the density function of Eq.(4) gives the probability that 
the paca. veccon velocity at any point in the gas will have a magnitude between R and 
R + dR : 


The same density function also represents the probability that a bomb will hit at a distance 
between R and R + dR from a given point if it is initially aimed at a point whose distance 
from the given point is P. The mean square aiming error is represented by Yo» the distri- 


bution being assumed Gaussian ‘*”) , 


156 


DEFINITION AND EFFECT OF DETECTION 


A detector is defined as any device whose instantaneous output is a function 
of the envelope of the input wave only. Thus 


R 
y= A) = F(v) (5) 


where yis the output of the detector and v is the normalized amplitude of the envelope. 


If P/V, is replaced by a, Eq. (4) may be written 


24a? 
dP = ve ? I, (av) dv, v>0 (6) 
dP = 0 v <0 

Eq.(5) solved for vis v = g(y) . : (7) 


If v is eliminated from (6) and (7), an equation of the form 
dP = f(a,y)dy (8) 


is obtained for the probability density for the normalized voltage at the output 
of the detector which has the characteristics given by Eq.(5). For example, if y 
= v?/2, then Eq.(8) becomes 


dP = e-%-* I, (2vzy)dy, y>O (9) 


dP = 0 y <0 


where a?/2 has been replaced by x. The quantity x may be identified with the power 
signal-to-noise ratio, commonly used in radar literature. 


EFFECT OF VIDEO AMPLIFIERS 


Since a complete radio receiver usually has one or more stages of video ampli- 
fication following the detector, it would seem that one would want to calculate 
the probability density function for signal-plus-noise at the output of the video 
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amplifier. This can be done theoretically, as has been shown in an excellent paper 
by Kac (#5), but the mathematical labor is great. On the other hand, it has been 
shown experimentally‘*®) that the signal threshold is practically unaffected by 
the video bandwidth until it becomes less than about % of the IF bandwidth. Since 
video bandwidths less than 4% cf the IF bandwidth are quite uncommon in practice, it 
appears best for the sake of simplicity to assume the video bandwidth infinite in 
all the work which follows. 


When the results have been computed, assuming an infinite bandwidth, it will 
be possible to modify them in an approximate manner so that they become valid for 
any video bandwidth. This is explained on page/213, under the title "Collapsing Loss". 


PROBABILITY OF DETECTION WITH NO INTEGRATION 


The calculations necessary to determine the probability of detection when 
exactly one correlation interval is available are quite simple compared with the 
case where the output over many correlation intervals is available, and hence the 
former case is taken up first. In a pulsed system this corresponds to using a Single 
pulse, while in a c-w system it is equivalent to observing the output for a time 
t = 1/Af, where Af is the over-all effective bandwidth. In either case this amounts 
to observing the receiver output for one correlation interval. If the output exceeds 
the bias level, the signal is observed or detected (see pages 9-14 of RA-15061, 
A Statistical Theory of Target Detection by Pulsed Radar ‘?®), hereafter referred to 
as No.1, for complete definitions of detection and bias level). 


It will now be shown that the probability of detecting a given signal x is 
independent of the detector finction, everything else being held constant and only 
one variate being taken from the density function of Eq.(8). The false alarm time 
has been defined as the time in which the probability is % that the noise alone 
will not exceed the bias level (Eq.(15), No.1), but it will be best here to keep 


things general and denote this probability as P,» rather than as 4%. Eq.(15), No.1, 
then becomes 


P= 1-Pi =P (10) 


where the subscript N denotes the number of variates and I’ is simply an abbrevia- 


tion. From Eq.(8), 
F(®) 
Pe AR f(0,y) dy (11) 


% 


where the symbol y, is now used for the bias number. Then the probability of detec- 


‘tion is 
F(@) 
ae f f(a,y) dy (12) 


% 
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but since y = F(v), or v = g(y), Eq.(11) may be written 


re) _v _8°(y5) 
r= i ve*dv=ee 7 : (13) 
& 


(¥,) 


and 
/ 1 
e(y,) = V2 log. - (14) 


Therefore Eq.(12) becomes . 


P = ve #7 I,(av) dv (15) 
V2 log, 1/7, ; 


which is independent of the detector function. 


The integral of Eq.(15) must be evaluated by approximate methods. This function 
will appear in several places subsequently, and is defined as* 


eee 
Q(a,8) = { ve §*  Ilav)dv . (16) 
p 


Footnote on Q Functions 


* It does not appear possible to express the Q function in terms of a finite number of 
known functions. The Q function is similar to Lommel’s functions and in fact can be ex- 
pressed as 

a® +5 


Qa.) =1-e * [iv,(-i6% ing) -u,(-i6%iap)] 


where U, and U, are Lommel’s functions of the first kind. This identity may be proven 
using the definite integrals given in Watson(1), pages 540 and 541, especially Eq.5 of 
page 541. By successive integration by parts, the Q function may be expanded in infinite 


series giving é 


2 +, r=0/a r 
Sie PE i I 
Q(a,A) = ¢ bal: (ap) 
or 
_ at +8 iz r 
Q(a,A)=1-e * =°(4) I (ag) . 
r=y 


The similarity of the first of these expansions to the series for U,(w,z) given in Eq. (1), 
page 537 of Watson, is interesting. A simple expression for Q(2,2) analogous to Eqs. (9) 
and (10), page 538 of Watson is 


Q(a,a) = Z [: see 1, (@*)| (Continued on next page.) 
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In terms of this notation, Eq.(15) may be written 


P, = Q(a,v2 log 1/T,) . (17) 


This is the probability of detection only if 7, = 1/f, where 7, is the time available 
for detection. In general the probability of detection is given by Eq.(18), No.1, 


ae Gel Gaal 0 hd (18) 


which follows from the definition of detection given on page 9, No.1. The double 
subscript notation P, fae is used here in place of P’. If N = y, Peg is written 
simply as P,. 


Ae (R/R, ) can be calculated by means of Eqs.(18), (17), the tables of Q, and 


the simple Cals ation 


R i 
Ry” V4 Va 


(19) 


(see Eqs.(10), (11) and (12), No.1). 


Footnote on Q Functions (Cont ’d) 


which is weer in special cases. An asymptotic expansion for Q which is of value is 
given by Rice (48), page 109: 


(4-a) 2 
Ae OL 1+(B- 
Q(a,f) = {1 = ¢71(p-a)] + i e ‘ car + late | 


2n 


where ¢ “4(T) is given by the error function of Eq. (100). This expression is most useful 
in the region where af >> ] and a >> | p-a|. 


The Q function is a special case of the incomplete Toronto function described in the 
footnote on page 182.| The relation is 


Q(a,4) =1-T, (1.02) . 
a v2 


The Q function is graphed in Figs.13 and 14. 


A table is available in Ref.47 but the intervals are too large to be of general use. Project 
RAND is computing an extensive table of the Q function which will be published as a sepa- 
rate report. 
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A very good approximation for the quantity I), which appears in Eqs.(10) and 
(17), may be derived from Eq.(10) by writing 


1/n' 1/n' log P, 1 1 
Pion ite gad ha log, p (20) 


which is valid when n‘ >>], a condition nearly always true in practice. Eq.(10) 
then becomes 


1 
log. B : (21) 


If FP, = 1/2, as is assumed in all the curves in No.1, 


0.693 
UG Brenan ual (22) 


Eq.(17) may consequently be written 


Eo 
" 
Lo) 

NI 


Ls 
BR) °¥4-60 log, m + 0.730) - (23) 


As an example, let R/R, = 0.595, andn = 10*. Then Pos Q(4, 4.37), which has the 
value 0.410 from the table given in Ref.47. Note that this is a point on the graph 
of Fig.l, page 22, No.1. 


GENERAL CASE — INTEGRATION OF N INDEPENDENT VARIATES 


If the output of the receiver (or filter) can be observed for a length of time 
much greater than one correlation period, it is of advantage to integrate the out- 
put. The simplest concept of an integrator is a device which linearly adds the 
voltage output of N samples from the detector. The time elapsing between samplings 
must be at least one correlation period, in order that the samples may be considered 
to be independent. If the sum of N variates* of signal-plus-noise exceeds the bias 
level calculated from the probability density function for N variates of noise 
alone, then the signal is said to be detected. 


* Readers with some statistical experience will recognize that here is a case of testing a 
statistical hypothesis. It is known that the n observations y,, yg, ---- y come from 
a universe whose density function f(y,a) depends on the unknown parameters a; it is required 
‘to decide, on the basis of these observations, which of the two values a, or ag is a better 
estimate for a. If a, is the true value of a, let p, be the probability of making the 
mistake of choosing a, as the correct value; similarly, if a, is the true value, let p, be 
the probability of choosing a,. Suppose p, = .05. Then a statistical decision method can 
be devised for which p, = .05 and for which P, will be less than for any other method with 
the same p,. See, for example, Kendall, vol.2, pp. 272-275 (6). 
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The integrator may take the sum of the squares of the N variates, or, in general, 
the sum of N variates where each variate has been processed by some general func- 
tion. As long as the same weight is applied to each variate, the integrator will be 
called linear. The function which the integrator applies to each variate will be 
called the law of the integrator. Any nonlinear integrator will be inferior in opera- 
tion to a linear integrator with the same law and would ordinarily never be used 
intentionally in practice. Cathode ray tubes are nonlinear, however, and thus fall 
short of other types of linear integrations. 


The law of the integrator acts in exactly the same way as the law of the detec- 
tor. Thus, if the detector output is y = (F(v) as given by Eq.(5), the integrator 
output is ’ 


N N 
Y = Ls = at aCade (24) 


It is obvious, as far as the theoretical problem is concerned, that the only 
function of importance is 


wv) = d[F(v)]. - . (25) 


There will be an infinite number of combinations of > and F functions which will 
produce the same function W and hence the same theoretical results. In all the work 
that follows, the output of the combination of integrator and detector for one 
independent variate. will be called y = ¥(v), and the sum of N variates will be 


¥ = Ys (26) 


1 


The symbolic solution for the case of N variates corresponding to Eq.(15) for one 
variate is not too difficult to obtain. The starting point is Eq.(8) for the proba- 


bility density function for one variate. The characteristic function for this 
distribution is 


c= f fa, yet dy . (27) 


@ 


The characteristic function -for the probability density function for the sum of N 
independent variates is then simply 


Cy = (G,)" (28) 
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and 


© 
avy ae 
dP, = av f Cy (a, w) e7 1Y (29) 


‘= 


or 
dP, = G(a,N,Y) dY . (30) 


Corresponding to Eq.(11) is 


Dace Ip G(0,N,Y) dy (31) 
, 
and to Eq.(12), 
P= I G(a,N,Y) d¥ . (32) 
Y, : 


If ¥, is eliminated from Eqs.(31) and (32), there results a solution for Py as a 
function itt, , N, and a, which is the desired result. 


It is found in most cases that the integrations required in Eqs.(27) to (32) 
are not possible in terms of known functions. 


THE SQUARE LAW DETECTOR WITH N VARIATES 


It seems, by a process of trial and error, that the best possible function 
for W(v) in Eq.(25) to produce integrable functions in Eqs.(27) to (32) is 


wv) = Av2 ay. (33) 


Though this represents a square law for the combined detector and integrator law, it 
is usual to think of it as representing a square law detector coupled with a linear 
law integrator. 
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In Eq. (33), the only effect of the constant A is to multiply the bias level 
Y, in Eqs.(31) and (32) by A. The value of P, in Eq. (32) is independent of A. It is 
Been ane to let A = 1/2, or y= v?/2. By direc substitution from Eqs.(6) and (27), 


“ @ 
C, = i e--* I, (2vzy )ePY dy (34) 
0 


where x = a7/2 and p = iw, 


This integral may be obtained from pair 655.1 of Campbell and Foster‘7). In all 
pairs taken from Campbell and Foster it is necessary to replace p by -p, since they 
use €?8 for the first integration. As long as the same notation is used in both 
directions, the order of signs is immaterial. In order to avoid confusion, the minus 
sign will be used in the exponent in the first transformation and the plus sign in 
the second transformation. Thus all of the characteristic functions which appear 
hereafter are really C(-p) rather than C(p). In this way there is direct agreement 
with the Campbell and Foster tables as well as with tables of the Laplace trans- 
form. Equation (34) becomes 


pa Ba 
C, = eae Mehta (35) 


The characteristic function for the sum of N variates is then simply 


‘ an es Basie. (36) 


By means of pair 650.0, Campbell and Foster, the probability density function is 


y 
dP, = (<=) eur (a/Nx¥)aY Y>0 (37). 
0 


Y<0 


Graphs of this function are shown in Figs.1-7. The density function for noise alone 
(x = 0) is found most easily from pair 431, Campbell and Foster, to be 


(38) 
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BIAS LEVEL FOR SQUARE LAW CASE 


The bias level, Y,, is by Eq. (31) 


@ 
\ y-1.-¥ 
7 L (N-1)! (39) 
b 
The incomplete gamma function, as defined by Pearson‘®), is 
uvpt1 
‘I(a,p) = if a ne (40) 
0 Px 
In terms of this function, Eq. (39) becomes 
Y, 
Paarl (,n-1) ; (41) 


The tables of I(u,p) extend to p = 50, and values of the function are given to seven 
places. Thus, for N < 50, and n’ < 10°, the bias level } may be obtained directly 
from Pearson’s tables. Other methods must be evolved for N > 50 or n’ > 10%. The 
normal approximation to Eq.(39) is unsatisfactory for N less than several thousand 
because of the fact that the integral is over a region which is far out on the tail 
of the curve. This can be seen from the Gram-Charlier series which will be taken 
up presently. 


The integral of Eq.(39) may be evaluated directly by successive integration 
by parts to give 


yo , 
oe Ad he BN xed hn CNELINE 2a : (42) 
Nn ~ (y-1)! b Y; 


In the regions of interest ¥,>N>>1. The series in the brackets may be approximated by 


Ope ee ee es x : (43) 
A | Nel 
Y 
b 
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in this region so that Eq.(42) becomes 


see -/, 
- aie ’ ; NY,e ; (44) 
v” N-1\ N!(¥,-N+1) 
CN-E aap 
b 
By the use of Sterling’s approximation for N!, Eq.(44) reduces to 
pe cement tos, 
exp|-¥,+N\1 +log, 7 
pe oe eT ee) 
een (Y, -N+1) 


Though the expression looks more cumbersome in this form, it is actually much simpler 
to use in calculations than is Eq.(44). Substituting for I) from Eq.(21) gives the 
expression 


login = 0.24 + 4 logy N + logy, (¥,-N+1) (46) 


Y 
+ 0.434 (¥,-N) - N log, + 


Graphs of this function are shown in Figs.8 and 9. For N = 1, the exact expression 


for ¥, from Eq.(39) is 
¥ = 2.3 log, n + 0.37 (47) 
whereas Eq. (46) for N = 1] reduces to 


¥ = 2.3log,,n+ 0.45 . (48) 


The difference is seen to be practically negligible. 
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CUMULATIVE DISTRIBUTION FOR N VARIATES OF 
SIGNAL PLUS NOISE — SQUARE LAW DETECTOR 


Knowing the required bias level for a given false alarm interval, it is now 
_ necessary to integrate the density function of Eq.(37) from this value to infinity 
to give the probability of detection for a signal of strength x, thus 


@ ied - Y-N% 
P, al ev N* 7 (2VNaY)aY . i (49) 


This integral is not soluble directly in terms of well-known functions. The order 
of the Bessel function* can be reduced in steps of 1 by successive integrations by 
parts, so that the last remaining integral is of the type given by the Q function 


* The following are some of the useful identities concerning the modified Bessel functions of the 
first kind: 


Ez) = A= -iP J (iz) = ue es AG [ gece’ OAT ea Aagemees ----| 
r=0 r! (ntr)! N! 2(2n+2) 2 °4(2nt2) (2n+4) 
2 2 
‘as r4 Zz 
I,(2) =1l1+ wa ur w 3 eeoe 


Asymptotic expansion: 


ge Ja BV i occas 
I (2) Poy aal [ + Jrhn® 4, Uir4n )(9-4n --] 
n V27%x 


1! (8x) 2! (8x) 


i (s* et=cosW (sin y)?" dy 


n 7” 
Saar aes 
vi T(nt5) 
zI/(z) = nl,(2) + 21,4 ,(2) = onl, (2) + 2I,_,(2) 
E,(3) = 1, (2) 

Sf 2"1,_ (dz = s"I (2) 

J 2,(2)d2 = 21,(2) 

2 = = 

2n 7 (2) = 1,_4(2) - Tyg, (2) 

Sf e*1,(x)dx = ze*{I,(2) - 1,(2)] 

J e"*I,(2dx = xe~*[I (x) + I, (2)] 

J e*I, (ads = e* (1-2) I(x) + sf, (2)] 
S e-*1, (ade = e°* (Cita) 1,(2) + 2, (2)] 


Relations between the J, functions and the hypergeometric functions will be found in the foot- 
note on p., 175. 
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of Eq.(16). An easier way to arrive at the same result is by the use of the charac- 
teristic function. To get the cumulative distribution from -© to Y of any density 
function, it is only necessary to find the anticharacteristic function of C/p, where 
C is the characteristic function of the given density function (see pair 210, Campbell 
and Foster). Thus from Eq. (36), 


bod e-Nz Nz PY; * 
ai ]i= ome @ Pe ld é (50) 
Mare f pips ge : 


-@ 


The term 1/p(pti)*" may be expanded in a series 


1 a 1 1 1 


po nil I Sehr is gh 2 A le Poa abaraa od aN ba Sb rT 
p(p+1)"% = p(p+i) (p+1)?— (p+)? (p+1)% ks) 


The mate of the first termof the series, by pairs 210, and 655.1, Campbell and Foster, is 
Y, 
om f e~ I, (/Nxy)dy : (52) 
0 
The first two terms of lia are thus 


Y 
- f * e-y-"e1 (a/Nzy) dy (53) 
0 


@ _v?+2Nx 
= f ve 2 I, ( w/2Nzx ) dv 
V3¥, 


= Q(VINx 2%) 


using the definition of.Q from Eq.(16). All the succeeding terms may be obtained 
by using pair 650.0 Campbell and Foster. 


Nx r-1 
1 


+ 2 
Mate of a is(X) e-YI (2VNxy) (54) 


— 


* As in Campbell and Foster, f is here used in place of w/2n or p/2ni. 
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From Eqs.(53) and (54), 


1 


rm 


Py = (ane, vam) +o ™ SRB)" oven) - (55) 


a 
Nx 


r=2 


This form of solution for By is practical for numerical calculation only where N is 
less than about 10. 


The characteristic function in Eq.(50) can be expanded in another manner using 


r=0 


1 1 


p(p+i)* De, (p+1)? Ve 


This leads to an expression for F, complementary to that of Eq.(55) of the form 


r-1 


~Y,-Nx Peery. 2S 
Prentice 62° tig oy (=) I, (2VNx¥,) . (57) 
r=Nt+1 Nx 


By equating (55) and (57), one obtains one of the known expansions for Q given in 
the footnote on page 159. Equations (55) and (57) may also be obtained directly from 
Eq.(49) by repeated integration by parts. Equation (57) may also be obtained di- 
rectly from Eq.(55) by means of the identity 


eh (2) (58) 


given in McRobert ‘*), page 32, and one of the known series for Q. 


For the special case ¥ = Nx, the function Q of Eq.(55) is simply 
1 see BY, 
Q=5 E e 1,(2%)| (59) 


(see footnote, page 159),' and Eq.(55) becomes 


1 -2Y,) 1 
En gt. © +133,(2%) +1, (2%) + 1,(2%) ---- Iy-1(2%)] : (60) 
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‘This formula is useful for checking special points for values of N around 10 or below. 


None of the methods developed above are suitable for calculating if for large 
values of N. ; 


In the next section the general method of Gram-Charlier series is developed, 
which will be useful in a number of succeeding problems concerning distribution 
functions over. large ranges of variation of N. 


EXPANSION OF FUNCTIONS IN GRAM-CHARLIER SERIES 


The function $(y) is defined by 


1 reat 
POD 7a=e ; (61) 
The Hermite polynomials may be defined by the relation 
i GU 
? (y) "me (y) (62) 


where the superscript 1 stands for the i** derivative with respect to y. The ¢ 
functions and the Hermite polynomials are biorthogonal, that is 


i; Hyd (y) = 8, = 0,84; (63) 


= 1,t=j 


Therefore it 1s possible to expand any reasonable function in a series of the form(®) 


fly) = Sa dity) . (64) 


t=0 


The coefficients a, may be evaluated in a manner analogous to the Fourier series 


methods by multiplying both sides of Eq. (64) by H(y) and integrating from -@ to ©, All 
terms drop out but one, giving 


Rs tay es 
a..aF =r if H; (y)fly)dy_. (65) 
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It is usual to make the substitution t = y - y/o before making the expansion, thus 
causing the second and third terms of the series to vanish. The notation is 


YORsu= the average value of y, or the first moment 


o? = at Nae the variance 


ves the n*® moment of the distribution. 


Vy, -f vrone (66) 
Equation (64) is replaced by 
Ey) a : padi (67) 
and Eq. (65) by 
c, fa coacnee ey ee rye : (68) 


It follows at once from Eq.(68) that c, = l/o, c, = c, = 0. The moments about the 
mean, or the central moments, are defined by 


KM; = f (y-y) fly) dy X69) 


foe) 


and the standard moments about the mean by 


Gee (70) 
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The coefficients c, in Eq. (68) may be easily written in terms of the a’s. The first 


’ few are 


C3 = Uehas 


iz) 
w 


1 
4 Tie 
SES a 
c st MSs 10a, ) 
1 
Cc = a! (a,-15a, +30) 
c, ® a (a,-21la,+ 1054, ) 


1 
c, = at (a,-28a, +2 10a,-315) 


1 
= pert (2,-36a,+37 8a,-1260a, ) 


(71a) 


(71b) 


(71c) 


(71d) 


(7le) 


(71f) 


(71g) 


Formulae for the y’s in terms of the v’s can be obtained directly from Eq. (69), giving 
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(72a) 


(72b) 


(72c) 


. 2 
He = Ye - 5UY,Y, + 10K, - 10v,v; + 47 (72d) 


- & ree 3 sag 6 
pe = Ve 6uY, + 15Y,r; 20v,y + 1Sv,w - SYY . (72e) 


Continuations of this series are obvious. 
The process of obtaining the Gram-Charlier expansion is now evident: 


. Find the moments of the distribution. 

. Obtain the central moments from Eq. (72). 

. Obtain the standard central moments from Eq.(70). 
. Obtain the coefficients from Eq. (71). 

. Write the series for f(y) from Eq. (67). 


Op won 


It turns out that the best grouping for the terms of the series of Eq. (67) is 
different from the natural sequence ‘®), Such a regrouped series is termed an "Edgeworth 
series" and is actually used in this work. The grouping used by Edgeworth is 


0 (73a) 
0.3 (73b) 
0, 3, 4, 6 (73c) 
0 ni 4 612 Si9des 9: se (73d) 


This means that if the 0 and 3 terms are used as the first approximation, the addi- 
tion of terms 4 and 6 gives the next order approximation, and so forth. 


MOMENTS OF SIGNAL PLUS NOISE, SQUARE LAW DETECTOR 


The moments of a distribution may be obtained by using the characteristic 
function as a moment generating function*. Thus 


: dt 
= (-))'(|— 74 
v, = (-1) tel) | (74) 


ane 


* Kendall, p. 54 (9) 4 
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In the case of the distribution function for the sum of N variates of signal plus 
noise with a square law detector, the characteristic function is given by Eq. (36), 
and the moments are 


i| @ eM epti 
Sit eee (75) 
i dp* (p+1)% p=0 


Though the first few moments may be obtained by direct differentiation, it is better 
in this case to expand in a McLaurin’s series and obtain the coefficient of p'/i!. Thus 


Nx 
+1 


ght Loge op labat Moats Esau Wa lieeee (76) 
(p+1)% = (p#1)% (p+ #2 (pi )at2- at | 


The coefficient of p'/i!is, by direct expansion of each term in Eq.(76), 


(N+i-1)! (N+i) | (N+i)(N+ti+1) (Nx)? 

j) pci B i sctitees ance 
Anni ND OND) n° | ae 
,(N+i-1)! ; 

=(-1) “wey! iF, (N+i,N,Nx) (78) 


where Ae is the confluent hypergeometric function.” Thus the moments are 


* The following are some of the useful relations concerning the confluent hypergeometric function: 


Tic) "5° _Tlatr) or +H) 
iF, (aee,2) = 5 utr) ara we: gre oe ' 
T'(a) reo r! T'(ctr) , c c(c F (a,b,c, ¢) 


Asymptotic expansion: 


2 
aA, (an es zy Leen: f + (ia) (c~a) za] 


(a) 2°~* z1! 
Kummer’s first transformation: 
fy (% c,z) = eek: (cma, ¢, ~z) 
Kummer’s second transformation: 


Pee 
ah, (a 24,22) = e of, (a+ =) 


Recursion relations: 


af (atl, ¢, z) + (a-c) if (a-l,c,z) = (z+2a-c) Ade (a, ¢, z) 


ac F (atl,c,z) + (c~a)z F (a, c+1,z) = c(atz) F { (a, ¢, 2) (Continued on next page.) 
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f (N+i-1)! 


Vv; ani! e7N* Fi (N+i,N,Nx) 


oF 


(N+i-1)! 
Vv. & 
(N-1)! 


F,(-i,N, +~Nx) , 


Hypergeometric function (Cont ’d) 
He (atl,c,z) + (1-c) ALA (a,e-l,z) = (atl-c) F (a,c, 2) 
c,F (a-l,c,z) + 2,F (a,ctl,z) = ¢,F (a,c, z) 
(a-c), F, (a-l,c,z) + (e-1) mA (a, c~1,2z) = (zta-1) fi (a, c, z) 


(c-a) 2, F (a, ctl, z) + e( e-1) Ff (a,c-l,z) = e(zte-1) Ff (a, c, z) 


4 if (ae. 2) a af (ath, ct1, z) 

Relations between hypergeometric functions and other functions: 
hy (Gree 2) = e* 
AA (a, at1,-z) = az™° ie e~*e2-1g¢ = 2° That IA, a-1) 


using Pearson's notation for the incomplete gamma function. 


= 77,2 ee 
F, (1, atl, z) = e7*z EOE eet a 1) 


F (3.4, -2*) = cs ot aa erf z 


ae 
af (on 1,2) = Li (2) 
(original Laguerre polynomial) 


ni T (atl) LI? (2) 


F (-n, ati, z) = 
11 (atl¢n) ” 


(generalized Laguerre polynomial) 


vil, 2+1,-2) = Teepe , (-z) 


arnt, (#) + 21, (9] 
F (4,1,-2) = =... ($) + 21, 1,(4)] 
n($2.-2) = ¢"#, (f)-4(8) 


' 
dol= 
= 
_ 
1 
x 
~— 
or 
wie 
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(79) 


(80) 


Eq.(80) being obtained from Eq.(79) by Kummer’s first transformation. The first 
four moments are 


v, = N(1+x) } (81a) 
v, = (Nx)? + 2Nx(N+1) + N(N+1) (81b) 
v, = (Nx)? + 3(Nx)?(N+2) + 3Nx(N+1)(N+2) + N(N4#1) (N42) (81c) 
v, = (Nx)* + 4(Nx)3(N+3) + 6(Nx)? (N+2)(N43) (81d) 


+ 4Nx(N+1)(N+2)(N+3) + N(N+1)(N+2) (N+3) 


The generalized Laguerre polynomial L{°\(z) is defined by 


T(atltn) 


(2) = —— 
bh fe n'lP(atl) 


olen, ati ele (82) 
Comparing (80) and (82), it is seen that the moments expressed in terms of the 
Laguerre polynomials are 


v= itLt(-Nx) (83) 


Another generating function for these polynomials is available through the relation 


& n 
LPG) a e722 - (e-2 2"*a) j (84) 


The moments about the mean may be expressed in terms of the moments about the origin 
by means of Eqs.(72a-e), resulting in: 


ie eed (85a) 


HW, = 0 (85b) 
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by = 2Nx + N = N(2x41) = o? (85c) 
pe, = 6Nx + 2N = 2N(3x+1) (85d) 
My, = 12(Nx)? + 12Nx(N+2) + 3N(N+2)  (85e) 
He = 120(Nx)? + 20Nx(5N+6) + 4N(SN+6) . (85f) 


A generating function for the central moments may be obtained by multiplying the 
generating function of Eq.(75) by e?”1 giving (see pair 207, Campbell and Foster) 


‘ = +n (Nx+tN) 
pb, = (-1)'e7%* at eres f (86) 
: dp* (p+1)% p=0 


The moments of Eqs.(85a-f) are most easily obtained by logarithmic differentiation 


in Eq. (86). 


The standard moments about the mean are obtained from Eq.(70), and are 


a, = 1 (87a) 

a, = 0 (87b) 

a, #1 (87c) 

ae 2 Eade Lunes (87d) 
NY2¢2x+1) 2 
6(4x+1) 

= oT ae 87 
ss ada * N(2x+1)? (87e) 
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“There is an approximate method of computing the significant part of a which is 


based on the fact that ¢, of Eq.(7ld) is always nearly equal to eZ (see page 
259, Fry‘®)). Thus : 


. 2 
a,x 15a, + 10a; - 30 (88) 
or 
10 (108%?+78x+13 ) 
a, ———— (87£) 
N(2x+1) 


For noise alone, the moments are given by 


;-7])! 
Vv. = Cibo es (89) 
ta (N= 1e 
and the central moments by 
(N+i-1)! ie ey 
K; = OE Mae a i * (90) 


Equation (90) was obtained from Eq. (86) by putting x =0 and expanding in a series. 


THE GRAM-CHARLIER SERIES FOR THE SQUARE LAW CASE 


The coefficients of the series may be obtained by use of Eqs.(7la-d) since 
the standard central moments are now known (Eqs.87a-f). They are: 


c¢-8 l (9la) 

ecm ec, #0 (91b) 

cfeper wa EL (91c) 
3N¥2 (2x41)? 
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c = 4x+1] 
4 4N(2x+1)? 


(3x+1)? 
CC ——————— 
6 18N(2x+1)°3 


From Eq.(67) the required series is 


Qa. 


epee 


where 


YY 
t= = ; 
foe N(1+x) 

o = VN(1+2x) 


and the c’s are given by Eqs. (9la-e). 


Ze, $8) + c, $9 (t) + c, $4 (t) + cg $8(t) -=== | 


(91d) 


(9le) 


(92) 


C93) 


(94) 


Note that the grouping .of terms is according to the Edgeworth scheme given in 
Eq.(73). Note further that as N tends to infinity, all the coefficients go to zero 


except c,. Thus 


as N >, In terms of N and x 


[y-n(1+x)]? 
fia ens CY Pee sc ae 
V27IN(1+2%) 
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(95) 


(96) 


Eq.(95) is precisely a statement of the central limit theorem, and the derivation 
given is essentially a loose proof of the theorem. 


The cumulative distribution is easily obtained from Eq.(92) by means of the 
simple relation 


fea = g(t) (97) 
giving 


Bi Hl fly)dy = ZU-$UT)] - PT) - G(T) - 6 A(T) ---- (98) 


N 
¥, 
where 
Y- vu 
b 
T-—— (99) 
and 
1 ie see 
s(T) =—— 7 da . 
f(T) al e a (100) 
aT 


The function ¢71(T) is tabulated in the W.P.A. tables‘®). This differs from the 
definition given for ¢-'(y) in Fry, page 456 but is used here because of theW.P.A. 
tables. 


The series of Eq.(98) was used to calculate all the curves of Figs.1-50, No.1, 
with the exception of the cases where N = 1]. In most cases the first two terms of 
the series are sufficient, though in some regions of small P four terms are needed. 


SAMPLE CALCULATION 


Assume N = 10, n = 109 
From Fig.8, or Eq.41, ¥, = 30.0 


R 
Let ils 1.0 so that x = 1.0 
0 
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From Eq. (94), 


Y= N(1+tx) = 20.0 


o = VN(14+2x) = 5.48 


‘From Eq. (99), 


30.0 -20.0 
ae EOE 1. 830 


$71(T) = £71(1.828) = 0.9325 


4(1-¢°T)) = 0.0338 


From Eq. (91c), c 


From Eq. (91d), 


From Eq. (9le), 


p? (1.828) 
$3(1.828) 
$5(1. 828) 
c,$7(T) = 
c, 93 (T) 
c¢,°(T) 


= 4X4] 
4 4N(2x+1)? 


(3x41)? 
c= SACLE ATES. 


18N(2x+1) 


wae B24) ni: 081 
3N¥7(2x41) 


2 
2 
= 0.0139 


= 0.0033 


0.174 


(See p. 218 for references on 
tables of the derivatives 
of the error function.) 


= -0.470 
= 0.990 
-0.0141 
-0.0007 


+0.0032 


P = 0.0338 + 0.0141 + 0.0007 - 0.0032 = 0.0452 


This point, P = 0.045, R/R, = ], may be found on Fig.20, No.1. 


INTEGRATION LOSS, 


SQUARE LAW DETECTOR 


It is of interest to express the effect of noncoherent integration as a loss 


with respect to coherent integration 
loss as 


u 


(43) | 


L. = 10 fg, 
1 


This may be done by defining the integration 


Nx, 
(100a) 
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where 


N = number of pulses integrated 


x, = required value of signal-to-noise ratio to produce given probability of 
detection for N = 1. 


x, = required value of signal-to-noise ratio to produce the same probability 
of detection for N = N. 


Thus L, is a function of P and n. However, it turns out that the dependence on P 
and mn is very small. 


In the case of coherent integration, x, is always equal to x,/N, so that L; = 
0. With noncoherent integration, x, is always greater than x,/N, so that noncoherent 
integration is never as efficient as coherent integration. The results of calcula- 
tions are given in Figs.10 and 11. One observes that the dependence of L, on P and 
n is quite small. Thus by means of the graph of Fig.12, which gives x as a function 
of P and n for N = 1, and any one of the curves of Fig.10, it is possible to obtain 
a fairly accurate value of x for any P, n and N. 


GENERAL CURVES OF THE CUMULATIVE DISTRIBUTION FUNCTION 


The integral of Eq.(49) is a function found in other applications than the one 
discussed in this paper. It is desirable to have graphs of this function available 
in general form rather than the specialized form of Figs.1-50, No.1. The integral 
is a special case of the incomplete Toronto function* described by Heatley ‘*®) and 
Fisher ‘'7), which is defined as 


B 
T,(m,n,r) = arate lt ret Tartide . (100b) 
0 


Using this notation, Eq.(49) for the cumulative distribution function becomes 


Py = 1 - Te (2N-1,N-1,VNe) (100c) 


* In normal correlation theory, the quantity | 


n-2 
= G1 7301152 
B -4(B7.°*) 
dp-4iq) ¢* 3, (aps)athns) 
a0a" *) 


is given by Fisher(17) as the limiting form, for large samples, of the frequency element 
of the quantity B? = n,R* where R® is the sample estimate of the multiple correlation 
coefficient of a random variable y with other variables x,, Ho xoe= » nM, 1s the size 
of the sample, and By = ngp2 where p is the population multiple correlation coefficient. 
The cumulative distribution is 


f= Tp (Gee 4 nit) 
en 


and can be obtained from the curves in Figs.13 to 32. 
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The function plotted in Figs. 13 to 32 is 
T(2N-1,N-1,V7q) (100d) 
and q may be found easily from these curves for any values of Y,. N and x. 


THE LINEAR DETECTOR — WN VARIATES 


The linear detector is usually more difficult to deal with than is the square 
law detector. The distribution function for one variate of signal-plus-noise is 


dP =ve ? I,(av)dv . (10) 


In attempting to find the distribution for the sum of N variates by the method 
of characteristic functions, the immediate trouble is that the characteristic function 
of Eq.(101) does not seem to be obtainable in closed form. To give an idea of the 
difficulty involved, the characteristic function for one variate of noise alone is 
obtained as follows: 


a” © v . 
C = ae jen: seem (102) 
0 


This is pair 903.3, Campbell and Foster, and may be evaluated directly by completing 
the square or by forming a differential equation, giving in either case 


3 
Cine) aie 2 eee (103) 
xy pe “erfe Te 


or in terms of w 


aw? a x? = 
C= ae ae” 4 { e2dx + iV (104) 


To raise this expression to the N** power and then obtain the anticharacteristic 
function is practically hopeless. 
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The distribution function for the sum of two variates of noise alone is ob- 
tainable by use of the convolution theorem, giving 


2 2 
agp Ha v . fy? 
2 > e (2, ca very dy (105) 


and the cumulative distribution is also obtainable, giving 


fo) 2 2 
2 deg ee 
nef fA vegeta a Teestag git ove (106) 
af 


However, when N > 2 there seems to be no closed solution corresponding to Eq.(105) 
or (106). Since these cases are for noise alone, the signal-plus-noise situation 
must be attacked by other means. 


It turns out that if the moments of the distribution for one variate are known, 
the moments of the distribution for the sum of N variates may be found directly. 
Formulae are given, for instance, in Cramer, page 345, (1°) for the first few central 
moments, which are 


PP = Nu, (107a) 
uy = Nu, (107b) 
uN = Nu, + 3N(N-1)p2 (107c) 
ee a 10N3 43 + 15N 33 + ISN? yu, - 45N? U5 : (107d) 


The corresponding coefficients in the Gram-Charlier series then become 


a 
3 

Pe: “31yiv72 (108a) 
a‘ ~3 
10a? 

ce Sine (108c) 
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The a’s in Eqs.(108a-c) are the central standard moments for one variate. Note 
that, in the square law case, if N is put equal to 1 in Eqs.(85c-f) and the resulting 
H's used in Eqs.(107a-d), the w’s for N variates are correctly given. If a moment 
generating function can be found for the case of N variates, then it is immaterial 
which method is used; but in the case in which such a function is not available, the 
Eqs. (107) must be used (or some method essentially equivalent). 


To handle the linear detector it is now sufficient to find the moments for one 
variate only. Rice'?®), page 107, gives the required expression as 


vas 2i/2r (143) F, (-4,1,-x) : : (109) 


Rice also gives the first two moments as 


; = (Ee? (+2012) + 21,(2) (110a) 


Vv, = .2(14+x) (110b) 


To calculate v one needs to know the function, #(-3/2,1,-%). This may be obtained 
by use of the recursion relation 


a,F,(at1,c,2z) + (a-c) F,(a-1,c,z) = (2at+z-c) F(a, c, z) (111) 


by putting a = -4, c = 1, z = -x. The result is 


x 
ue 
V, = 2v, (2+x) pine =1,(=) (110c) 
also 


V, = 4(2+4x+x7) (110d) 
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The corresponding central moments are 


pe ** o? = Q(1l+x) - vu (112a) 


x 


av? - 2v, (1422) Ze ?1,(5) (112b) 


Ms 


My = 4(2+4x+x7) - 34 - 4v, [(a-ady, - /Eer*1,(2)| < “Cii2e) 


The standard central moments, and then the c’s of Eqs.(108a-c), are directly obtain- 
able from these formulae, though the process is somewhat tedious due to the cumbersome 
form of Eqs.(112a-c). The functions v4 to 4% are shown graphically as a function of 
x in Fig. 34. 


To obtain the bias level ¥ for the linear detector for N > 2, one can use the 
G.C. series for noise alone. Setting x = 0 and v= V772 in Eqs.(112a-c) gives 


Hy = O72 = 2-5 = 0,429 (113a) 
ML, = \Ein-s) = 0.1772 (113b) 
2 
wy, = 8 - 2 = 0.598 (113c) 
and 
oes = = 0.632 (114a) 
Oo 
0.1053 
c, vv? (114b) 
and 
be 
a, ae eS. 26 (115a) 
= 
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3 0.0108 


: ; (115) 
0.00555 
Di eceraeeras tat aia 

N (115c) 
The cumulative distribution function is now equated to es giving 
0.693N 1 -1 021093 0.0108 0.00555 
Cees 5 [1-4 (T)] fae 0) ms Sa = apr Tae ---(116) 
where 


For any given nm and N, Y,may be found from Eq.(116) by trial and error methods. 
If an approximate value of T is found by neglecting all but the first term in Eq. (116), 
a more accurate value obtained by Newton’s method is 


f (1) 
Neh LGR) 


CHA) 


It is better, however, to plot Fq.(116) giving n as a function of T and N from 
which is finally obtained the bias level graph of Fig. 35 showing ¥, as a function 
of n and N for the linear détector. 


Since for finding the bias level it is necessary to know the distribution 


functions only for large values of the argument, it is possible to find an ap- 
proximate solution valid in this region. Consider a distribution function given by 


dP = ve * dv (117a) 


for v going from -ato+o. The N™ convolution of this function will be nearly 
the same for large values as if (117a) went only from 0 to ®, because the large 
values in the sum of N variates are most probably produced by addition of large 
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values of every variate, and for large values (in fact for all positive values) 
the two distribution functions are identical. The characteristic function of Eq. 
(117a) is given by pair 710.1 of Campbell and Foster to be 


2 


Pe 
C= -/27 pe’ . (117b) 
For the sum of N variates 
No 4, cup? 
Coe Ca 1) 29)? pee =, (117c) 


The probability density function is obtained from pair 740.2 of Campbell and Foster as 


x. 2 
GMa cry 
FT ns ee Dy (sz) y>>1 (1174) 
N 2 


where D, is the parabolic cylinder function of order N. In terms of the derivative 
of the error integral as defined in Eq. (62), 


N 
(277)? y 
OR Re re $"(Te)dy Pee en a (117e) 
N 2 
Note that for N = 2, Eq.(117e) becomes 
2 
aL 2 
Ree es alte 


Referring to Eq.(105), the exact expression for this case, it is seen that Eq. (117f) 
can be obtained by neglecting the first term and replacing erf y/2 by 1, both of 
these approximations being very good if y >> l. 


The approximate cumulative distribution is easily obtained from Eq.(117e) by 
direct integration and gives 


N 


Py, (22) av-1 (1) (1178) 
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The bias level is easily obtainable from this expression by equating it to [, and 
solving for } by means of the tables of ', or by plotting graphs. The method is 


not very practical for N > 20 since suitable tables do not exist. 


It is interesting to note that no such approximation as Eq.(117g) is obtainable 
for the square law case. 


Graphs of the probability density functions for signal-plus-noise have been ob- 
tained by numerical convolution for some selected cases and are shown in Figs.36 to4l. 


RESULTS OF THE LINEAR DETECTOR CALCULATIONS 


The difference in results for the linear and square law detectors turns out 
to be so small that extreme accuracy must be used in the calculations to show the 
relation in its true form. One such comparison graph was calculated and is shown 
in Fig.42. Also, in Fig. 43 is shown the difference in db in the two cases at P = 
0.50. The two are identical at N = 1, the linear law becomes better by a maximum 
of 0.11 db at N = 10, the two are again equal at N = 70, and the square law then 
becomes better and asymptotically exceeds the linear law by 0.19 db as N ~© having 
reached 0.16 db at N = 1000. These results show conclusively that there is little 
to choose between the linear law and square law as far as theoretical signal threshold 
1s concerned. 


EXPANSIONS IN LAGUERRE SERIES 


In certain cases, particularly for low values of N, the Gram-Charlier series 
may not be the best-suited type of expansion for distribution functions which are 
zero for all negative values of the amplitude. For low values of N, a suitable 
expansion for such functions is the following: 


fly) = 5 a, e-%y7 L3G) (118) 


1=0 


where L*(y) is the generalized Laguerre polynomial defined by Eq.(82), or by 


Zy-@ dt ; 
L2(z) = = epee lies (119) 


The orthogonality relation which makes the expansion possible is 


- 4 Le 9) 
rein f en? ALi (z)Le(z) = 8; (120) 
0 
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(see Copson‘*), page 269). 


Thus, from Eqs.(118) and (120), the coefficients are determined by 


a BY 
Sp a : (121) 
< ea Ee eee i 


Note that Eqs.(118), (120) and (121) are analogous to Eqs.(64), (63) and (65), re- 


spectively, for the Gram-Charlier expansion. 


Let. a new variable t = y/G. Then 


t 
8 


fly) = g(t) = ¢,e"? toh F(t) (122) 


t=0 


where 


aE ceo a eee ae a(t et ee 
“i l(atit1) if bite eet os reeeee) Af Le(Z)f iy) 2 (123) 


The first few Laguerre polynomials are 


Li (zjee ad (124a) 
LS (2) = lta-z (124b) 
2L9(z) = (a+1) (a2) - 22 (at2) +z? (124c) 


6LG5(z) = (at+1)(at2)(a+3) - 3z(a+2)(at+3) + 32? (a+3) -23 . (124d) 


There fore 


= ea ae 1+ we 
AiavRrE Rta hes oe (125) 
Lares Bc) Let 
c, atari) (ar2) B (a+2) =] (126) 
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Since there are two arbitrary constants, @ and £, in the expansion of Eq.(122), it 
is possible to make c, = ¢, = 0 by a proper choice of a and £. These relations are 


easily determined by equating Eqs. (125) and (126) to zero and solving simultaneously. 
The results are 


1 "1 
a =——>-1-=—>-1 (127) 
Dead o 
v-v? 2 
pee - > | (128) 
1 1 
and 
1 cl 
Bp ae ot SS ee (129) 
& BI(atl) gar(2.) : 
o2 
c= C= 0 (130) 
Paley el 97 pha (131) 
= ———|/— (a -_—= 5 
“3 “ Br(a+4) |B? B3 


The coefficients past c, are so complicated that the whole value of this type of 
series seems to depend on the fact that the first term alone is often a good ap- 
proximation. This approximation is 


Y, ee Wy af 
| DBS ae GS Bry. dy eeeG132) 
and the corresponding cumulative distribution function is 
4 y y? 
p= f fly) dysesdce: Le Al (133) 


o? 
y 
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where I is the incomplete gamma function as defined by Eq. (40). 


There is a striking analogy between Eq.(132) and the corresponding normal 
approximation. In both cases the distribution for the sum of N variates is simply 
obtained by multiplying both y and ao? by N. As N +, both the normal approximation 
and Eq.(132) approach the true distribution (and each other). In any particular 
case, however, the convergence properties of one approximation will be more useful 
than the other. 


In the square law case, for x = 0, 4, = N and o? = N, Substitution of these 
values in Eq.(132) gives 


dP = : e-YyN-ldy . (134) 


Note that this is the same as Eq.(38), the exact expression. Thus in this particular 
case the first term gives the whole correct result. The third coefficient from 
Eq.(131) is easily shown to be zero, as all the following coefficients will be. 


In the square law case where x # 0, y = N(1+x) and o? = N(1+2x). Substitution 
of these values in Eq.(132) gives 


N(1+x) 2 
Ate ee eee 
ah RI Se ial ei)» (12) (135) 
N(i+x)? 1+2x 
(1422) Sy ay 
and from Eq.(133), 
y N(1+x)? 
Pa Se eS ‘ 
cross 1+2x ] eee 


A comparison of the particular case N = 3, x = 1 is shown in Fig.44. Curves are 
given for the exact distribution function (Eq.(37)) and the two approximations given 


by Eqs.(96) and (135). 


For the linear case with x = 0, 4 = W7/2 and o? = N(2-7/2), the cumulative 
distribution is, from Eq. (133), 


peal hee rien (137) 


YN(2-77/2.) ° 4/m-1 
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OTHER SERIES APPROXIMATIONS 


It is theoretically possible to develop still other series approximations for 
the various distribution functions. For, instance, it might be thought advantageous 
to use a sum of terms of thetype y*e"-” /2 particularly in the linear case. While 
this turns out to be possible, even the first coefficient is so difficult to calcu- 
late that the process is impractical. 


METHODS OF INTEGRATION INVOLVING 
SUBTRACTION OF NOISE 


Certain practical difficulties arise in maintaining the bias level at the 
correct value in an electronic detector, particularly if the number of pulses inte- 
grated is large. The trouble may arise from fluctuations in amplifier gain, the 
bias supply, or the noise level itself. 


A solution of this problem is to have the gain of the amplifier, or the bias 
level, or both, controlled by some sort of average value of the noise output. Ob- 
viously the time constant of the control device must be neither too long nor too 
short. One scheme which has been used is to subtract a pulse known to consist of 
noise only from each possible signal-plus-noise pulse* (see paragraph 3, page 11, 
No.1). Thus in the absence of a signal, the average value of any number of composite ~ 
pulses will always be zero, and the required bias level will be comparatively low. 


DISTRIBUTION FUNCTIONS FOR COMPOSITE PULSES 
OF SIGNAL-PLUS-NOISE MINUS NOISE 


When a noise pulse is subtracted from each signal-plus-noise pulse, the theo- 
retical distribution functions will be entirely different from previous cases. The 
square law case is the only one that can be treated in any reasonable fashion. The 
distribution function for one variate of signal plus noise is given by 


dP = e** 7 (2/xY) aY (138) 


and the characteristic function is 


+1 
a eeceen ala (139) 


* This subtraction can be accomplished by means of a gate which operates at double the 
repetition frequency. On every other gate only a noise pulse of reversed phase goes through 
the integrator. 
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Subtracting a positive noise variate is equivalent to adding a negative noise variate. 
The distribution function for a negative noise variate is 


fy ger Bena (140) 
3 = 0 Y>0 
and 
Gee (141) 
1-p 


To obtain the characteristic function for the sum of a variate from the distributions 
of Eqs.(138) and (140) it is only necessary to take the product of the characteristic 
functions given by Eqs.(139) and (141), giving 


e7= ee 


e Pt (142) 


1-p? 


This is the characteristic function for one so-called composite pulse. The characteris- 
tic function for the sum of N composite pulses is simply 


e7 Nz ae , 
= —————- eP 143) 
Co a oa? 
In the case of noise alone (x = 0), 
: (144) 
(1-p2)% 
and the anticharacteristic function is, by pair 569 Campbell and Foster, 
1 
dP ane 
1 en ie Y\d 
N  Vvr(N-1)! A Red : ses 


where K,_, is a modified Bessel function of the second kind and is given by the 
finite series 


r=N-1 
ee (N+r-1)! 
Ky.y(2) = Vote? Ss ST (146) 
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The cumulative distribution for the sum of N composite noise variates may be found 
by use of the series (146) and term by term integration. However, for N greater 
than 3 or 4 the process rapidly becomes impractical. 


Again, it is necessary to find moments and proceed by means of Gram-Charlier 
series. For noise alone, the moments are easily found from Eq.(144) to be 


De set peas todd (147) 


(wei-1)4(i)! 
Ue pe Se ae FNS eh he TEER 


(N-1) 1(3\ 


in particular, 


LL, = 2N = o? (148) 
2, = 12N(N+1) (149) 
Hd, = 120N(N+1)(N+2) (150) 
and 
a, = 0 (151) 
guheus ie (152) 
30 
Pe ee (153) 


The only coefficients different from zero in the first six are ¢, and ¢, to the 


4 
order of 1/N. 


c= = . (154) 
8N 
‘Thus 
psa (=) i ‘(s) as (155) 
ah ale Jan) * an® \Van 
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and similar to Eq.(98) is the cumulative distribution 


P, = ho (A) - <3] “--- (156) 


The bias number is found by setting this expression equal to Ty and plotting ¥, 
as a function of n and N. Results are given in Fig. 45. In the special case N = 1, 
the cumulative distribution function is simply e~'’? for Y > 0, and the bias number 
is obtained from this expression rather than from Eq.(156). The anticharacteristic 
function of the general case, Eq.(143), may be obtained by use of the convolution 
theorem, pair 202, Campbell and Foster. Let 


Cie 
pt eos (157) 
1 = (p+1)% 
1 
a) nee 158 
a (ep wee 
Then from Eq. (37), 
N-1 
2 
G, = eo eNsT  (2/Nxy) y>0o (159) 
Nx # 
= 0 y <0 
and by pair 525.2, Campbell and Foster, 
Gye" 0 y>o (160) 
a ee 
(N-1)! — 
Applying the convolution theorem gives 
N-1 
e t-Nx sod y 2 
aye a | (ge) OTDM ey, (QVNay dy ¥>0 (161) 
‘ Y 


For Y < 0, the lower limit of the integral in Eq.(161) is 0 rather than Y. 
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PROBABILITY DENSITY FUNCTIONS FOR Y < 0 


To evaluate the integral in Eq.(161) when the lower limit is zero is straight- 


forward but tedious. First one evaluates the integral 


N=1 


Geta Ce eaees VNxy 
f = : (<2) e7*Y I, (2Y/Nxy) dy 


(162) 


by use of characteristic functions in a manner entirely similar to that used in 
Eqs. (74) to (80). The characteristic function of the function of Eq.(162), with 


k= 0; is 
Nx 
er 
(p42) 
and 


Nx 
e 7(N+k-1)! Nx 
ICR) = ayrowek 4 Webs) 


Then by expanding (y-Y)%~!, one obtains the coefficient of 


(N-1)! (-y)e be 


fp See 
k! (N-1-k)! 


and from Eqs.(164), (165), and (161), 


@ Enna pies wes 
k! (N-1-k)1 


dP, = 
| (N-1)! k=0 


or 


2 ae rs (- _Nx 
3 k=n-1(N+k-1)! Ff k,N, sheets Y <0 
(N-k-1) $k! Qt 


dP, = dY 
x (N=) ees 
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(163) 


(164) 


(165) 


(166) 


(167) 


In terms of Laguerre polynomials, using Eq. (82), 


l= k= 
(eke) 1a GYR Po” Yes Ore, (168) 


k=0 


The first few polynomials are given inEqs.(124a-d). Some special cases of Eq.(168) are 


ae 

Neat) fi idPysdie Y <0 (169a) 
etx 

N=2 dP, = ayS—(1s2-Y) Y<o (169b) 
4 2 
1-3 


9 2 
seaxe2Z (2022) + v2] y-<0 (169c) 


The cumulative distributions for Y < 0 may easily be obtained by integrating (169a-c). 


Obviously, the expressions in Eqs. (167) and (168) are practically useful only for 
small values of N. 


PROBABILITY DENSITY 
FUNCTIONS FOR Y > 0 


To find a general expression for Eq.(161) giving the distribution function 
when y > 0 is a task of tremendous proportions. Consider, for instance, the special 
case N = ]. Equation (161) becomes 


© 
dP, = ayers e 7) I,(2vxy )dy 120-43 (170) 
: 


By means of the substitution y = v?/4, this becomes 


Plas BY pe ae Se 
dP, = ra e 7? I,(wx)dv y >On, (17h) 
avy 
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This can be expressed in terms of the Q function defined by Eq. (16). 


y-% 


z 
aP.*s" a¥ =O Veav¥).)  ¥ 30.8 (172) 
Eq.(169a) was 
"5 
Seer TY 0 (169a) 


Thus, for the case N = ], the whole distribution function is described by Eqs. (169a) 
and (172). A graph of this function for various values of x is shown in Fig.46. 


Note that if x = 0, Q(0,2VY) = e7?’, and Eq.(172) for Y > 0 reduces to 


ev! 
OP ieee end Lact 22240 (173) 
and from Eq. (169a) 
et 
a Boe ee (174) 


when x = 0. Thus over the whole range of Y 


-lr| 
oP dy (175) 


which checks Eq.(145) when N = 1. 
For N = 2, Eq.(161) becomes 


ee 
dP, = aver { (a) athe” I,(2/txy)dy  Y>0. (176) 
Y x 


This integral may also be evaluated in terms of the Q function. The process requires 
a large number of integrations by parts and is very time-consuming. The result 
turns out to be 


e Y= 2x 
4 


dP, = dY S (142-1) (az, 2/7) + }¥5(o/ ; Lymm)! 


Ya. (177) 
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‘ 


This equation is already so complicated as to be nearly useless. Thus it was not 


thought worth while to seek a general expression of this type for arbitrary N when 
K>:.6, 


Note: If x = 0 in Eq.(77), it reduces to 
en! 
dP, = dyY<-(1e¥) —-¥>0 (178) 
which may also be obtained from Eq.(169b) by substituting -Y for Y. 


CUMULATIVE DISTRIBUTION FUNCTIONS 


The effort in Eqs.(161) to (178) has been concerned with obtaining the proba- 
bility density function for N variates of signal-plus-noise minus noise. To find 
the cumulative distribution functions exactly is difficult, especially for Y positive. 


A case which can be solved, however, is that for N = 1]. For Y negative the 
answer is simply obtained from Eq.(169a) and is 


ee 
e 
Por lais 5 ) i Va (178a) 
For Y positive, using the result of Eq.(172), 

ae 

read (ig 
P= =A eY Q( Vx, 27y) dy ¥30-,. (178b) 

Y 


ae 
Since the value of A at Y = 0 is, from Eq.(178a), 1 - e 7/2, Eq.(178b) may be re- 


written as 
eee Gad 2 5, hare e%Q(Vx,2/7y) dy (178c) 


but from the definition of Q in Eq.(16), 


@ v2 +x x © 
Q(vx,2/y) = if we 8 I, ( wx) de = 2e a be I,( 2vxz) dz , (178d) 
y 
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Replacing Q by its defining integral in Eq.(178c) gives 


x 
2 . i 
Pee 2s ste etf dyer f e-?? J (2vxz) dz ‘ 
0 of 


Integration by parts is now used, letting 


fo) 
oe f e-?? T (2vxz)dz 
y 


dv = eYdy 


du = -e° I (2/xy) dy 


v = ev 
x 
pales ala. oe? (aay e 
UME tye 2e ‘4 o( 27xy)dz - > 
ya se 
= Q(x , 2/Y) ie oe 
Thus 
e? J 
P=l1- 7 #1 27¥) + of vdu 
0 
or 


y-% 


2 Y 
Piom ies —— Ave ,2V¥) - sf eY I,(2vxy ) dy : 
0 


(178e) 


(178f) 


(178g) 


(178h) 


(1781) 


(1783) 


The integral term in Eq.(178j) is just 1 - Q(v 2x,/2Y), and the final result is 


yo 


P, = Qv2x,/2Y) - £— vz, 277) Y>0 


2 


B 
For x = 0, Q(0,8) = e?, and 


Y 
ee ewe ee) ree 
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(178k) 


(1781) 


agreeing, as it should, with the result obtained from Eq.(145) by letting N = 1 
and integrating. For Y = 0, Qa,0) = 1, and 


A a ar (178m) 


agreeing with Eq.(178a) when Y = 0. 


The bias number for use with Eq.(178k) is obtained by 


-y, 
0.693 e° 
i Cerra as oe (178n) 
or 
¥ = 2.30 log, n- 0.327. (1780) 


In Fig. 47 is shown a graph comparing Eq. (178k) with Eq. (23) for n = 10°, where 
P is plotted as a function of x. 


Though it might be possible to calculate the cumulative distributions for N 
> 1 by a method similar to that used for N = 1, it would be very tedious. Therefore 
resort is made to Gram-Charlier series, as before. The moments are directly obtainable 
from the characteristic function given in Eq. (143), 


dt “Nx rial 

. e e 

. « (-])! | 179 
v. = ( Was fice | ae (179) 


There seems to be no readily obtainable expression for v,in closed form. The first 
s1x moments obtained directly from Eq.(179) are: 


v, = Nx (180a) 
v, = (Nx)? + 2Nx + 2N (180b) 
v, = (Nx)? + 6(Nx)? + 6Nx(N+1) (180c) 
v, = (Nx)* + 12(Nx)? + 12(Nx)?(N+3) + 24(N41) + 12N(N+1) (180d) 
vz = (Nx)? + 20(Nx)* + 20(Nx)9(N+6) + 120(Nx)?(N42) + 60Nx(N+1)(N+2) . 


(180e) 


202 


UV, = (Nx) + 30(Nx)> + 30(Nx)*(N+10) + 120(Nx)>(3N+10) 
+ 180(Nx)?(N+1)(N+6) + 360Nx(N+1) (N+2) 


+ 120N(N+1)(N+2) . (180f) 


‘The corresponding central moments are: 


HM, = QNx + QN= QWN(i+x) =o (181a) 
é, = 6Nx (181b) 
M, = 12(Nx)? + 24Nx(N+1) + 12N(N+1) (181c) 

H, = 120(Nx)? + 360(Nx)?(N+3) + 360Nx(N+1)(N+2) + 120N(N+1)(N+2) - 
(181d) 

The central standard moments are: 
3 
a, =—— > (182a) 
V2N(1+x)? 
3(1+2x) 
29 182b 
2 N(1+x)? hee 
45 (3x7+3x+1) 
a, wis + Nien (182c) 
The coefficients of the series are, from Eq.(71): 
6 = - ——— (183a) 
2/2N(1+x)? 
1+2x 

= ——— 183b 
c, ENC ipa)2” (183b) 
e id dan Deal (183c) 


“6 ~ T6N(1+x)3 


Y 


The Gram-Charlier series for the probability density function is given by Eq.(93) where 


i. —— oe Nano ee eae) (184) 


and the cumulative distribution is given by Eqs.(98) and (99). Figures 53 and 54, 
No.1, showing the comparison between the ordinary case and the composite case, were 
computed using the Gram-Charlier series developed above. There appears to be no 
significant difference in the probabilities of detection for N between 1 and 10. For 
N between 100 and 1,000, the composite case gives an effective signal-to-noise ratio 
about 1 db lower than the ordinary case. 


ANOTHER APPROACH TO THE DETECTION 
CRITERIA-PROBABILITY THAT SIGNAL- 
PLUS-NOISE EXCEEDS NOISE ALONE 


The method of setting a bias level and calling any signal-plus-noise or noise 
alone which exceeds this level a signal is not the only possible way of defining 
detection. Another method is based on asking what is the probability that any given 
signal will be larger than any noise pulse during a given interval of time‘?!), The 
interval of time taken would logically be the false alarm time, as defined previously. 
In this time there will be n/N = n' independent groups of noise pulses. If the proba- 
bility that a single integrated group of signal-plus-noise pulses exceeds a single 
group of noise pulses is called A(x,N), then the probability that the group of 
signal-plus-noise pulses exceeds all of the n’ groups of noise pulses is simply 


P = [P\(x,N)]” (185) 


This probability is a little difficult to interpret properly. It means that if 
during the false alarm time a signal of strength x appears, it will have this proba- 
bility of being larger than any noise pulse group appearing during the same time. The 
difficulty is how to pick out the largest signal over a period of time, and what 
to do when many signals are present. These are reasons why the earlier detection 
criteria are thought to be superior, since they provide clear answers for the above 
questions. The criteria presented above may be of special value, however, when a 
target is known to be present. Such is the case when a target is being automatically 
tracked, and one wishes to calculate the probability that it will be subsequently 
lost due to the noise exceeding the signal. 


The probability density function for N signal- plus- -noise pulses minus N noise 
pulses has been indicated in Eq.(161). 


To obtain the probability that the sum of N signal-plus-noise pulses will be 


greater than N noise pulses it is only necessary to integrate Eq.(161) from 0 to 
©, It will be easier to obtain the probability that N noise pulses exceed N signal- 


204 


plus-noise pulses, however, since this requires the integral from -® to 0, and an 
expression is available for Y < 0 in Eq.(167). Thus 


_Nx Ny 
e? 0 ken-1 (N+k-1)!,F,(-k,N,- Na) 
ea d Y SES see Att We eka MY N-k-1 (186) 
ee one m » (Nok-1)tktowe "8? 


Now one substitutes z for -Y and interchanges the summation and integration signs, 
obtaining 


_Nx _Nx 
oF bem |(Wek-1)!,F, (-k OW, af Sete | 
0 


oe Ghia 1s) Une (N-k-1)!R!Qntk €187) 
The integral is simply (N-k-1)!, and therefore 
7 N-1 
Sa ire (N+k-1)! Nx 
ity 2 = ats Diane (188) 
Or in terms of Laguerre polynomials, using Eq. (82), 
“, 
e 2 R=N-1 ; Nex 
2 =nek pn-a ( _N% ; 
Ler Ds gra-k prt (-— ) (189) 


From Eq.(188) a more convenient form may be obtained by introducing a dummy index 
t and interchanging summation signs, leading eventually to 


Nx, Nx 5 
Se at ("F j eS _(N+k-1)!_ 
a2 2 $t(es- il (N+i-1)! (ki) PgNeet 


1=0 


(190) 


The outside summation in Eq.(190) is obviously a polynomial in x of the N-lth degree 
and with N terms. It is rather curious to note that if one puts t = 0'in Eq.(190), 
the following identity results: 


k=N-1 
(N+k-1)! 
atom ede RA EE: 191) 
: > (N-1)!R!2 
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In other words, the constant term in the polynomial is always unity. 


The first few cases for low values of N are: 


yucé 
A = 2° 25 (192a) 
R= =e (1+7 (192b) 
HO: i) 
1 el 9x 9x" 
Rese 15 (192c) 
eis ( 20x x7 a 
A = 3 e7** (1 aon tae ae (192d) 


(192e) 


l a 325x 575x 4375x° 625x4 ) 
- \+—— + ——_ + ———_ + 
256 1024 6144 . 32,768 


Obviously for N very large, these expressions rapidly become useless, and it is 
necessary to use the Gram-Charlier series of Eqs.(184) and (98). The lower limit 
Y is replaced by zero, giving for the series 


P= —[i-g(T)] + h2(T) - G(T) - ogG8(T)---- (193) 


ole 


where 


Y, N 
T crcl Fae (194) 


and c,, ¢, and c, are given by Eqs.(183a-c). A graph of P as a function of x and N 


is shown in Fig. 48. For very small values of P, more terms may be necessary in the 


series of Eq.(193). 


USE OF CUMULANTS IN OBTAINING 
-GRAM-CHARLIER SERIES COEFFICIENTS 


It is often much simpler to obtain the cumulants for a given distribution 
function rather than the various moments. The cumulants may be defined by 


di 
K, = Cvs log, c) (195) 
P an 
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where C is the characteristic function of the given probability density function 
(see pages 61-65 of Kendall ‘®)). The cumulants, except the first, are invariant 
with respect to a change of origin. Also, for the distribution of the sum of N 
variates, it is only necessary to multiply every cumulant by N, as is evident from 
the defining Eq.(195). The coefficients of the Gram-Charlier series in terms of 


cumulants are given on page 149 of Kendall. The cumlants in standard measure may 
be defined as 


K. 
U 


K; ica (196) 
o 


In terms of standard cumulants, the coefficients of the series are: 


Cf ag EN BOM SL (197a) 
K, 

cee, (197b) 
K, 

Cc, = a (197c) 
1 2 (197d) 

c = (Ket 10K5) 


The first term in Eq. (197d), K,, is omitted in the 0,3,4,6 approximation. 


Consider the square law case where, from Eq. (35), 


-% <a 
ro orcs (198) 
p+l 
x 
log, Cx#-«x +I - ln (p+tl1). (199) 
From Eq. (195), 
K, = (i-1)! (ix) Ree oe enon (200) 
For N variates, 
K, = N(i-1)!(ix+l) . (201) 
and 
ie ee peed Ee) (202) 


t 


o* aol : 
NA (2x+1)? 
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In particular, 


< 2(3x+1) 


= (203a) 
2 


N1/2(9x+1) 


6(4x+1) 


* * N( 2x41)? fee 


and it is at once evident that Coa es and c. obtained from Eqs.(197b-d) are iden- 
tical to the values given by Eqs.(92c-e) by means of a much longer process. 


In the case of a composite pulse of signal-plus-noise minus noise, the charac- 
teristic function is given by Eq.(142) and 


LOgeeG, #7 <4 € = Leg (1-p?) . (204) 


Again by means of Eq.(195) it is easy to derive, for N variates, 


K, = N(i-1)! [ix+2] i even 
= N(i-1)! (ix) iodd, 71 (205) 
or 
K, = N(i-1)! fizere(-1)4] | (206) 
and 


K; (i-1)! [ixei1+(-1)' 
Ko = ; a]. : (207) 
o at Lt 
N2 [2(x+1)]? 
Special cases are: 
6x 


OG ise sme a RT (208a) 
Ni/2 [2(x+1)] 


r]w 


3(2x+1) 


K, = NCwe1)? (208b) 


and again by Eqs.(197b-d) the coefficients are seen to be the same as given by 
Eqs. (183a-c). 
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In a case such as the linear one where the characteristic function cannot be 
obtained, the cumulants are still useful and may be found from the moments v; by 
means of the formulae at the bottom of page 63 of Kendall. The first few are: 


Kk =Yy (209a) 
Ky = Vy = (209b) 
Kose Va tty, oleu: : (209c) 
Ky = Uy - hu, vy 8y a4 love = br... (209d) 


The K,; are now obtained by multiplying by N and dividing by o', The coefficients 
are then obtained as before by Eqs. (197a-d). 


BEST POSSIBLE DETECTOR LAW 


It is of considerable importance to know whether there may be some detector 
law which will give results which are appreciably better than the linear or square 
law cases which have already been considered. 


The problem may be stated as follows: 


These are available N samples 


dP, = ve 2 dy (210) 
or the distribution 
_yPta? 
dP, = ve 2 I, (av) dv (211) 


the former being the distribution of the envelope of noise alone, and the latter 
the distribution of signal-plus-noise. 


The probability that all of the variates v,--+--v, came from the second dis- 


N 
tribution is simply 


dP), = 4P,(v,)dP,(v,) ---- dP, (vy) (212) 


209° 


‘ 


whereas the probability that they all came from the first distribution is 


ar = dP, (v,) aP, (% ) woe dP, (vy ) (213) 


The ratio of dP, to dh, is the best measure of the likelihood that all the variants 
came from the signal-plus-noise distribution. It can be shown that any monotonic 
function of this ratio gives an equally good significance test. One arbitrarily 
picks a constant which the ratio must exceed to say that it shows that the variants 
came from the signal-plus-noise distribution. This constant determines the false 
alarm time. 


Taking the ratio of Eq.(213) to Eq.(212) and substituting values from Eqs. (210) 
and (211) gives 


we SE Raat (214) 


or 


a” i=N 


e ? |] I,(av,) 22 (215) 
t=1 


where \ is the constant which determines the false alarm time. 


Taking the log of both sides of Eq.(215) gives 


u=N 2 
> log, I,(av,) 2 log, A + i (216) 


i=0 

Note that nothing has been said in the foregoing discussion about integration. 
Now, however, Eq.(216) says that the best thing to do is take the log of I, of 
each variate, add these functions for each variate, and require the sum to exceed 


a certain value. Clearly this calls for a detector and integrator which has the 
combined law 


y = log I,(av) é (217) 


The meaning of this result is really quite remarkable (at least to one who is not 
a statistician). It says, in effect, that by having the sum only of N variates 
which have been subjected to the law y = log I, (av), one has as much useful in- 
formation as if the individual values of each of the variates were known (as far 
as determining to which distribution the variates belong)*. 


* If the two distribution functions to be distinguished are normal, then the simple sum of 
the N variates, or the mean, is the best criterion. In other words, a linear law would 
be the best if the envelopes of noise and signal-plus-noise were normally distributed. 
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Suppose that the signal strength is very small (which would make N large for 
any reasonable probability of detection). Then ‘I, (av) © 1 + a?v?/4 and 


ay? a? y? 
y = log I,(av) % log \1+ Pa cead (218) 


In this case, the square law is seen to be the best possible choice. If, on the 
other hand, the signal strength is large, I,(av) x e*’//27av and 


1 
y = log I,(av) = log Fea HT log 27av % av (219) 


27Tav 


Thus, for large signals (usually small N) the linear law is best. 


It should be pointed out that the results for the two extreme cases, square 
and linear law, are not very different (see Fig.42), and in practice a linear de- 
tector would usually be preferred on account of its relative immunity to saturation 
by large signals. 


In the case of a human operator it is difficult to say what law is used in 
the process of integration. Thus if a linear detector were used in the receiver, it 
is conceivable that the operator might mentally take the sum of the squares in his 
integration process, with a net over-all square law effect. 


SIGNAL-PLUS-NOISE MINUS NOISE — LINEAR LAW 


This case is of special interest because of the method which must be used in 
obtaining the solution. Since the characteristic function for the linear case cannot 
be found, it is necessary to determine the moments for a composite variate directly 
from the moments for the signal-plus-noise distribution and those for the noise 
distribution alone. 


Using a double subscript notation, in which the first index represents the 
number of the distribution function and the second index represents the order of 
the moment, the following formulae can be derived at once by successive differen- 
tiations of the product of the characteristic functions of the individual distribution 
functions: 


1 Gatos » lignes U (220a) 
My = Vig + Mag + 2M 1% 9 (220b) 
Vz = yz + yy + 3(%~M%_ + MQM) (220c) 
Vee se feu eA Ys Fey) 65 My (220d) 
Ve = Ye t+ Ug + 6(YysM%y + MYysYy) + 15(MQ%y + MyQ%y) + 20%3%, (220e) 
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The first set of moments are those for one variate of signal plus noise given 
by Eqs.(109) and (110a-d). The second set of moments are for one negative variate 
of noise alone. These are simply obtained from the first set of moments by putting 
x = 0 and multiplying the odd moments by -1. 


The details will not be given, since the results bear the same relation in 
general to the square law case as they do when a noise variate is not subtracted 
from each signal-plus-noise variate. 


USE OF SO-CALLED DETECTION CRITERIA 


Lawson and Uhlenbeck have made use of a quantity which is the shift in average 
value of a distribution of signal-plus-noise from that of noise alone divided by 
the standard deviation of noise alone, which they call the detection criterion. In 


symbolic form 


i 1S+N “1N ( 
= oy : 221) 


This quantity is also called the deflection criterion, and it is implied that it 
must be of the order of unity or greater to have a reasonable probability of detection. 


For the square law detector, using the results of Eqs.(8la-b) and (85c), the 
criterion becomes 


k = x/N (222) 


and for the linear detector 


xVN 


En 


k= 


= 0.957xVN (223) 


assuming x to be small. 


The object of these criteria is to show the variation in necessary signal- 
to-noise Tatio as a function of the number of pulses integrated. The results for 
k in Eqs.(222) and (223) may be derived. rigorously from the basic distribution 
equations if the central limit theorem is assumed to hold and for probability of 
detection equal to 0.50. 


However, it is found from the actual results presented in No.1, Figs.1-50, 
that the square root of N law given by the detection criteria is not closely fol- 
lowed, even for N as large as 1,000. If a law of the form 


k = xN? (224) 
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is assumed, the exponent 0 may be obtained from the data of Figs.1-50, No.1. The 
results are given in Figs.55 and 56, No.l. It is seen that 0 goes from 1.0 at N = 
1 to around 0,75 at N = 1,000. As pointed out earlier (page 182), 9 = 1 for co- 
herent integration. 


It has been said that the N“ law seems to fit observed data fairly well. It 
is the belief of the author that this is a coincidence that arises from the fact 
that the losses due to nonlinear integration by cathode ray tubes, and human operator 
losses, tend to just about equal the difference between N° and N%, so that the N% 
law actually seems to fit the observed data. 


It 1s rather interesting to note that if the detector law is assumed to be of 
the form y = uv", the detection criterion turns out to be 


nF!) 
pn avVN (225) 


avn! (21)? 


A graph of this function shows a very broad maximum of ] at n = 2. Thus this is a 
special case, showing that for large N the square law is the best of the particular 
class of functions v", This is not as general as the proof on page 210 which shows 
that the square law is the best of all possible functions for small x. 


COLLAPSING LOSS — INTEGRATION OF GREATER NUMBER 
OF NOISE VARIATES THAN OF SIGNAL-PLUS-NOISE VARIATES 


In many radar applications, an additional number of noise variates are integrated 
along with a given number of signal-plus-noise variates. Such is the case when 
three-dimensional data are compressed onto a two-dimensional presentation, or with 
a C scope where range is not shown. The loss so occasioned is called a collapsing 
loss (*®), An effect of the same kind is caused if the spot of a cathode ray tube 
indicator moves less than its diameter in a pulse length‘!!)., Again, if the video 
bandwidth is narrow compared with the IF bandwidth, the same sort of thing happens. 
All three effects are handled by assuming a given collapsing ratio, ~, which is 


defined by 


M+N 
N 


(226) 


p= 


where 


N = number of signal-plus-noise variates integrated 


M = number of effective additional noise variates integrated. 


In the case of loss caused by low writing speed of the cathode ray beam, the effective 
collapsing ratio is given approximately by 


d+sT 
POSE VS ee (227) 
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where 


d = spot diameter 


a 
i] 


writing speed 


T = pulse length. 


Where the loss is caused by a video amplifier, the equivalent defining equation is 


Fagin late op (228) 


where 


Bis = IF bandwidth (or total combined RF and IF bandwidth where RF 


amplification is used) 


B, = video bandwidth. 

Mathematically, the treatment necessary to take account of M extra noise variates 
is rather simple. It is only necessary to multiply the characteristic function for 
N signal-plus-noise variates by the characteristic function for M noise-alone variates. 
In the square law case, this results in 


no(4) 
oe oP EG) pti 


N (ps1 (p41) NP 


(229) 


It is apparent, by comparison with Eq.(36), that the results obtained for p = 1 
can be used directly to obtain results for any p. 


Care must be taken in obtaining the bias level, however. Without the M extra 
noise variates, the relation n’ = n/N is used to find the required signal-to-noise 
ratio, x. With the added noise variates, the number of groups of pulses integrated 
may or may not remain the same. In the case of video mixing, where the output of 
two independent radars is superimposed on the same indicator, the number of groups 
of pulses integrated is constant, which means that n’ is constant. 


In the other cases where the loss is caused by narrow video amplifiers, col- 
lapsing of coordinates, or slow writing speed, the number of independent groups 


of pulses integrated is reduced by the factor P.¢¢ so that n remains constant as is 
easily seen from the equations 


n= n'N (no loss) (230) 


n = (pn')(M+N) = n'N (with loss). (231) 


The collapsing loss is defined as 


c 


L, = 10 log, > (232) 
1 


where x, is the required signal-to-noise ratio with M extra noise variates, and 
x, is the signal-to-noise ratio required with no extra noise variates, such that 
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the probability of detection is the same in both cases. This fixed probability 
level will usually be taken as 0.90. 


The procedure, after finding x,, 1s to get the required bias from either Fig. 

8 or 9, depending on whether 7 or n’ is held constant, using PN as the number of 

variates. From the cumulative distribution functions graphed in Figs. 13 to 32, the 

value of x, is found by multiplying the finding x for PN variates to give P = 0.90 

and multiplying this value of x by p. The reason for multiplying by p is apparent 
on referring to Eq. (229). 


The results of the calculation are shown in Figs. 49 to 52 where L is plotted 
as a function of N for P = 0.90 and N = 10°. Also given are curves of @ defined by 


x 
eScrar Heke (233) 
x 


1 


It has commonly been said that 0. should be 1/26)» (*5), This statement is sometimes 
derived from the detection criterion given on page 212. 


From Fig. 57 it is seen that if n' is constant, 9 does approach 1/2 as N~> ©, 
However, 9 is much smaller for reasonably small N. In the case of n constant, the 
square root law is not even approached as an asymptote. 


It was found that the values of L and 9 are only slightly dependent on the 
original values of n and P. 


ANTENNA BEAM SHAPE LOSS 


It has so far been assumed that the antenna pattern was flat over the half-power 
beamwidth and zero elsewhere. In any practical case the beam shape may usually be 
approximated by a Gaussian curve which will hold fairly well out to + the beamwidth 
from the point of maximum gain. In the case of a searchlighting antenna, the re- — 
turned pulses will all fall at the same place in the beam, and if this does not 
happen to fall at the maximum of the beam, the loss may easily be taken into account 
by modifying the expression for gain used in Eq.(9), No.1 for calculating R, such that 


62 6? 
nano +22) 
G=G e is eee (234) 
where 


= azimuth angle between target and antenna axis 
= elevation angle between target and antenna axis 


half-power azimuth beamwidth 


a a ad ad: 
W 


= half-power elevation beamwidth 


If the antenna is scanning, the problem is entirely changed because the suc- 
cessive returned pulses will be of different magnitude. It is obvious that as the 
antenna scans past a target, pulses should be integrated out to some point where 
the principle of diminishing returns sets in. It is not too difficult to determine 
this point and to calculate the loss occasioned due to the beam shape as compared 
with the ideal case‘*®). A complete’ treatment which covers the general case of delay 
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of the received pulse relative to the transmitted pulse, off axis in elevation 
while scanning in azimuth, and random orientation of the pulse pattern relative 
to the antenna pattern is quite involved. However, the solution of some special 
cases has shown the general character to be expected of the results. 


The integration of pulses should be carried to about 1.1] times the half-power 
beamwidth. This figure is practically independent of the signal strength (range) 
and the number of pulses per half-power beamwidth. When the optimum number of pulses 
are integrated there will be an average loss over the ideal case which assumes 
constant gain between the half-power points. This loss is in the neighborhood of 
1.5 db and does not depend much on signal strength or number of pulses per hal f- 
power beamwidth. Since this loss is so small it was not considered worth while to 
reproduce all the detailed calculations here. 


It should be mentioned that special care is necessary when one considers rates 
of antenna scanning so fast that about only ] hit per beamwidth is obtained. In 
this case it may be expedient to make the receiving antenna lag the transmitting 
antenna to compensate for the time of travel of the pulse, or to step-scan, that 
is, move the antenna in discrete steps rather than continuously. 


In order to calculate the probability of detection in any case where the suc- 
cessive returned pulses have different signal strengths, it is necessary to obtain 
the over-all characteristic function by multiplying the characteristic functions 
for each pulse. Using this method it is not difficult to work out the needed results 
in any particular case. 


LIMITING LOSS 


If limiting occurs anywhere in the receiver, the probability of detection 
will be lowered, everything else being held constant. The video amplifier is the 
first place where limiting will probably occur. Let the limiting ratio be defined 
as the ratio of the limit level to the R.M.S. noise level. Limiting can then be 
represented mathematically by replacing the probability density function at the 
detector output by an equivalent function below the limit level, and a delta func- 
tion at the limit level having an area equal to all of the area of the original 
function to the right of the limit level. The moments can be calculated for these 
new functions (noise alone and signal-plus-noise), and the probability of detection 
found by use of the Gram-Charlier series as usual. The calculations are quite tedious 
and will not be reproduced here. The main conclusions are that if the number of 
pulses integrated is large, the limiting loss is only a fraction of a db if the 
limiting ratio is as large as 2 or 3, but if only one or two pulses are integrated 
the limiting ratio must be in the neighborhood of 10 to prevent a serious loss. 


Limiting in the output of the integrator can also cause a loss, but this loss 
is small compared to the loss caused by limiting of the individual pulses in most 
practical cases. 


EFFECT OF SIGNAL INJECTION ON PROBABILITY OF DETECTION 


It has been proposed that the minimum detectable signal can be decreased by 
the injection of an RF or IF carrier voltage that adds linearly to the received 
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echo and the receiver noise‘®!), The theory is that: the total signal will then be 
large compared with the noise, and thus the so-called modulation suppression that 
occurs in the process of detection with small signals will be eliminated. 


In such a process, the coherence of the injected signal with the received 
echo must be taken into account. If the target is moving, then the successive re- 
ceived pulses may be considered to be random in phase, so that the injected signal 
will necessarily be noncoherent with the echo. Analysis has shown that in this case 
the probability of detection decreases continuously as the magnitude of the injected 
signal increases, assuming a linear or square law detector. However, it can be 
shown that the best possible detector law starts to change radically as soon as 
the injected signal strength becomes comparable to noise. The analysis of proba- 
bility of detection when the detector function is altered to take into account 
the injected signal has not been completed. Preliminary estimates indicate that 
there will be only a small decrease in sensitivity in this case. 


It might be imagined that coherence could be obtained in a system using only 
one hit per target but having, say, 20 separate receiver channels with 20 separate 
injection oscillators having phases spaced 12 degrees apart. Thus, the return echo 
would be nearly coherent with some one of the channels. Theoretically, the improvement 
in this channel would be about 1] db. However, even this improvement would be just 
offset by the increased false alarm number due to the multiple channels, so that 
the over-all system improvement would be nil. It seems that there is no way to 
increase system sensitivity to moving targets by signal injection. 


There is some possibility of increasing sensitivity for stationary targets 
by coherent signal injection, but it is difficult to imagine a practical situation 
where such a method would be of any use. 


PROBABILITY OF DETECTION WITH MOVING TARGET INDICATION SYSTEMS 


The analysis of the probability of detection for MTI systems is quite compli- 
cated. It depends on the type of receiver (lin-log limiting or IAGC), the type 
of detector, and the characteristics of the storage device used. For a nonfluc- 
tuating clutter and no scanning noise, the effect of the clutter with or without 
the addition of a coherent oscillator is much the same as that of the injected 
carrier discussed in the previous section. If a suitable detection system is used, the 
sensitivity may be reduced by a small amount, due to the addition of the coho, perhaps 


by 1 to 3 db. 


The sensitivity of an MTI system for high probabilities of detection is further 
reduced due to the fact that the target may be moving at a speed differing from 
one of the so-called optimum speeds. This effect is quite complicated and is similar 
to that caused by a random variation of the cross section of a target with aspect. 
A method of quantitatively treating these problems has been developed and will be 
presented in detail in a future report. 


If there is a fluctuation component in the clutter, due either to the movement 
of the clutter itself or to the scanning of the antenna, the effect will be to 
increase the amount of noise at the receiver input. This can be taken into account 
by an appropriate adjustment in the value of the noise figure of the receiver that 
will change R, by the correct amount. 


217 


‘ 


TABLES OF THE DERIVATIVES OF THE ERROR FUNCTION 


In order to make efficient use of Gram-Charlier series, it 1s necessary to 
have a good table of the derivatives of the error integral (the ¢ functions of Eq. 
62). No satisfactory table was in existence at the time this report was written. 
Typical of the available tables ‘®) were 


Fry (°) n= 1(1)6, x = 0(.1)4 5 decimals 
Jorgensen nm = 1(1)6, x = 0(.01)4 7 decimals 


and an unpublished table of the W.P.A., giving 
n= 1(1)14, x = 0(.1)8.4 20 decimals 
RAND therefore decided to calculate a suitable table with the aid of its IBM 
equipment. This has resulted in a table of Hermite polynomials, as well as in the 
derivatives of the error integral, giving 


nvl=(TC1310, -~ =0€01) 12.0 6 significant figures 


A limited number of these tables are available at the present time. (RAND Document 
D-350, A Table of Hermite Polynomials and the Derivatives of the Error Function. ) 
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SUMMARY 


This report considers the probability of detection of a target by a 
pulsed search radar, when the target has a fluctuating cross section. 

Formulas for detection probability are derived, and curves of 
detection probability vs. range are given, for four different target 
fluctuation models, 

The investigation shows that, for these fluctuation models, the 
probability of detection for a fluctuating target is less than that for a 
non-fluctuating target if the range is sufficiently short, and is greater 
if the range is sufficiently long. | 

“The emount by which the fluctuating and non-fluctuating cases. differ 
depends on the rapidity of fluctuation and on the statistical distribution 
of the fluctuations. Figure 18, p.307, shows a comparison between the 


non-fluctuating case and the four fluctuating cases considered. | 
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SYMBOLS 


scale factor, Case 5 

scale factor, Case l 

incomplete gamma function 

number of hits integrated 

false alarm number 

probability of detection 

range 

range for which average input signal-to-noise ratio equals unity 
input signal-to-noise power ratio for a single pulse 

average of x over all target fluctuations 


normalized threshold 
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I. INTRODUCTION 


The probability of detection of a target by a pulsed search radar has 


been treated in considerable detail by J. I. Marcum? 2) 


for the case in 
which the amplitude of the returned signal pulses is not fluctuating. 

The purpose of this paper is to extend some of Marcum's results - mainly 
the computation of probability of detection* vs. range curves - to several 
kinds of fluctuating targets. (No claim to originality is made. Several 
of the equations given below have been derived by Marcum in hitherto 
unpublished work; some of these equations have also appeared elsewhere 
in'the literature. This report tig for the purpose of making the results 
more readily available than heretofore.) A general familiarity, on the 
part of the reader, with Marcum's paper will be assumed. Marcum's 
notation will be used throughout. 

Four different models of target fluctuation will be considered. The 
four models chosen for consideration are felt to be typical of situations 
which are likely to be met in practice, or, at least, to bracket a wide | 
range of practical cases. 

In applying probability of detection computations to actual cases, 
one should first attempt to analyze the Pietuntions of the actual target 
under consideration, and then choose whichever model (including Marcum's 
non-fluctuating model) appears to most closely approximate the actual 
target fluctuations. Or, one could consider the Ania target to be 
intermediate between two of the theoretical models. 

* Use is made of the term iprobabiiity of detection’ in order to conform 
to Marcum's terminology. Actually 'blip-scan ratio" is a more common 


term to use for the quantity which is computed. (The curves presented 
here all correspond to the cases which in Marcum's paper are labeled 


¥=N.) 
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One of the uses of the ensuing results is, that by comparing the 
results for different models, one can make some judgment as to the errors 
introduced by choosing the wrong fluctuation model. 

The four fluctuation models considered are as follows: 

CASE 1 

The returned signal power per pulse is assumed to be constant for the 
time on target during ea single scan, but to fluctuate independently from 
scan to scan. (This ignores factors such'as beam shape effect.) Expressed 
in statistical terms, the normalized autocorrelation function of target 
cross section is gicacianeens one for the time fn which the beam is on 
target during a single scan, and is approximately zero for e time as long 
as the interval between scans. This type of fluctuation will henceforth 
be referred to as scan-to-scan fluctuation. 

The fluctuations of target cross section are evidenced as fluctuations 
of signal to noise ratio in the receiver. The probability density for the 


input signal-to-noise power ratio is assumed to be: 


iM 7 AX/X * 
w(x,x) =—— e (x 20) (ra 
x 
where x = input signal-to-noise power ratio 
xX = average of x over all target fluctuations 
CASE 2 


The fluctuations are independent from pulse to pulse. This type of 
fluctuation will be referred to as pulse-to-pulse fluctuation. 


The probability density function is given by (I.1) 


* This formula also represents the probability density for target cross. 
section D> if x is replaced by 3s and x by L!. 
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CASE 3 . 
Scan-to-scan fluctuation according to the probability density 


be _“2x/X > 
[4 


w(x,x) = ear e (x 0) (1.2) 
x 


CASE 4 


Pulse-to-pulse fluctuation according to (1.2). 


It would be well at this point to indicate in which actual situations 
the various models would be likely to apply. 

As to the choice of probability density function for the fluctuations: 

Theoretically, for a target which can be represented as several 
independently fluctuating reflectors of approximately equal echoing area, 
the density function should be close to exponential, even if the number of 
reflectors is as small as four or five. Thus one would expect objects 
which are large compared with wavelength (and which are not shaped too much 
like a sphere) to fluctuate approximately according to the exponential 
density (I.1). 

On the other hand, a target which can be represented as one large 
reflector together with other small reflectors, or as one large reflector 
subject to fairly small changes in orientation, would be expected to 
behave more ‘like (1.2). 

The non-fluctuating model would apply to spherical or nearly spherical 
objects (e.g. balloons, meteors) at fairly large wavelengths. 

Most available observational data on aircraft targets indicates 

reement with the exponential density (1.1). More definite statements as 
to actual targets for which (I.2) or the non-fluctuating model apply must 


await further experimental data. 
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As to the choice between scan-to-scan and pulse-to-pulse fluctuation: 

The scan-to-scan model would apply to targets such as jet aircraft or 
missiles, for radars having fairly high pulse repetition rate and scan rate. 

Pulse-to-pulse fluctuation would apply to propeller-driven craft if the 
propellers contribute a large portion of the echoing area; or to targets for 
which very small changes in orientation would mean large changes in echoing 
avea, such as long, thin objects at high frequency; or to targets viewed by 
a radar with sufficiently low repetition rate. 

Most actual targets would probably be intermediate between the various 
cases considered. 

A comparison between Cases 1, 2, 5, 4, and the non-fguctuating case is 
given in Fig. 18 (for typical false alarm time and number of hits integrated). 

In all cases, it is assumed that there are on each scan N hits; after 
passage through @ square law second detector, the resulting N pulses are 
added and required to exceed a threshold 4% in order for detection of a 
target to occur.* The second detector output is, for mathematical convenience, 
assumed to be normalized as follows: detector output equals input envelope 
squared divided by twice the mean square input noise voltage. 

Formulas for probability of detection P, as a. function of x are given 
for each case in Section II. Curves (corresponding to Marcum's) for Pp vs 
range are given in Figs. 1-18. (Marcum's full range of parameters is not 
duplicated, but the means for doing so are given by the formulas in 
Section II.) Since the formulas in Cases 1 and 3 lend themselves to very 
convenient approximations, these cases are further discussed in dactten it. 
The derivations of the formulas in Section II are given in Section III. 

* @he actual beam shape is in effect being approximated by a beam having 
uniform gain over a finite sector, and zero gain outside this sector. 
In principle it is possible to take account exactly of beam shape in 
computing probability of detection. This is not done here, however, 
since it is thought that the effect of the aforementioned approximation 


is small provided the effective number N of hits per scan and effective 
xX are properly chosen. 
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II. FORMULAS FOR Py 


CASE 1 
The exact formula is 


Nee iis Py = eo | | 
Rani PUP iim 1 lee 
N-1 
ieee [=a 
(+ 3) i eescs =e exp | Ee (IZ.1) 
(3) 


where I is the incomplete gamma function. 


In most cases, this can be closely approximated by 


Tha GES “4 
Py = (1 + iz ) exp | | (II.2) 
For N= 1, (II.2) is exact. 

For N >1, in most cases of interest, Nx is several times greater 
than one. If this is true, and if the false alarm probability is sufficiently 
small, then, for N >1, the gamma function factors are very nearly unity. 
It turns out that in most cases of interest, the values of Nx and false 
alarm probability are such that (II.2) can be regarded as practically 
exact. (One must keep in mind, of course, that the assumptions Nx >1 and 
negligible false alarm probability mean that one cannot use (II.2) for 
indefinitely small values of Py: In most cases of interest, the applica- 
‘pility of (1.2) extends at least down to P, = one per cent. ) 


Now, for Nx > 1 the expansion of InP, as given by (II-2) is 


InP, = - os (X, -N+1) wars (x, ” aad 2 Feees (x, - Nh) toes (11.3) 
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‘oop : R 
This is @ convenient way in which to compute P,. Curves of P, vs p~ based 
: ° 


on (II.2) are given in Figs. 1 - 3. 
_ A good approximation to P, can be obtained by using only the first 


D 
term of this expansion, giving 
Py = exp =h{a,8) = : (Ir.4) 
x 
where 
Y7N+1 
h(n,N) = —¥ (11.5) 


Here n is the false alarm number; - is a function of both N and n.* 


a | 
Since = = (F-) » this leads to the result 
° 


Py = @ 
(II.6) 
R, g(n,N ) 
and | 
¥, “Nel ne et 
g(n,N) = N (Ti.7) 
In other words, the Py vs S curves are all » to @ good approximation, 


simply exp (aay with ~ = atete factor 1/e(n,N). Curves of exp (=u) 
and g(n,N) are given in Figs. 4 and 5. 

To see approximately the error introduced in going from (II.2) to 
(11.6): for each Py let Ry be the range computed by using just the first 
- term of the expansion (II.3); let Ry be the range computed by using the 


first two terms. Then**. 


* Curves of vs N and n are to be found in Ref. 2. In most cases, n is 
approximately equal to the false alarm. crag aaxided by the mune width. 
** - See Arpendie A for. derivation. 
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R, | : 
os 21 - E(n,N) np, (12.8) 

where N-l 
e(a,n) <2 5 (12.9) 
4(Y, -N+2 ) | 


Curves of €(n,N) are given in Fig. 6. 

Tais is valid for -16 €(n,N) 1n P, less than about .5. Referring to 
the curves of €(n,N) in Fig. 6, this means for Py greater than .Ol ~ .10, 
depending on the values of n and N. In almost all such cases, the first 
two terms of (II.3) are the only significant ones, so that (II.8) can be 
regarded as a correction between (II.6) and (II.2). 


The curve P., = exp (-u") is not necessarily the best curve to use in 


D 
(II.6) to give the best numerical results in approximating (II.2). It is 
nae to give this curve an average correction so as to obtain better 
results over the interesting range of the parameters n and N. 

Such a corrected curve, to be used with the scale factor g(n,N) to 
approximate (II.2), 4s given in Fig. 4. ‘The use of this curve with g(n,N} 
gives agreement with (II.2) to within about five per cent in R. tor | 


R 
6< 12 


° 
Py os Ol, 10 en = 10 » endl $n $1000. the agreement gets better as 


Pp gets larger, and is practically exact for Py 2 i. 
CASE 2 
The exact formula for Py is 
P,=1-I pad wer iy Eh (II.10) 
» (14x) JN 
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Reference 3 contains tables of the incomplete gamma function enabling 
the computation of (II.10) for N - 1 $ 50. Beyond this point an Edgeworth 


series can be used to compute Pp: 


pe k [ap en] - c,p(@¢xy - op ny - cg) (aa) 


where 
2 
s(t) = pe om || 
r 2 
gi(t) = To / exp (=) at 
. -T 
gi?) gs), gi?) are derivatives of ¢ (11.12) 
and 
¥, > N(1+x) 
fh (14x) 
0 = eI C= tes Ce = sae (11.13) 


R 2 
Curves of Py vs Ro for this case are given in Figs. 7 = 9. 


CASE 3 
The exact formula in this case is so cumbersome as not to be worth 
writing out here. (The exact density function is given p Eq. (III.18). ) 
The formula obtained by assuming the negligibility of false alarm 
probability, and S* >1, is 
Fd g 
Py = (1 + 25) ‘ Pe “fx 0-2) os il (12.14) 


N 
1+ l+> 


This turns out to be the exact formula for N = 1 and N = 2, Curvea 


of P, vs e based on (II.14) are given in Figs. 10 ~ 102, 
) 
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A close approximation to (II.14) is given by 


-N+2 +N+2 . 
Py = |2 +§ e = ) exp | -¢ (* . (11.15) 


1 R 4 
end since = = (F} » this gives 


R 
fe) 
Poa + ou") exp | - zu" } (12.16) 
R_ u 
Ro * <¢(a,N) 
where 
¥, “N42 1/4 

f(n,N) a yi (II.17) 


The scale factor f(n,N) is almost indentical to the factor g(n,N) of 
Case 1, the maximum difference being about 2 per cent for small N, and 
negligible for moderately large N. Curves of (2 + out ) exp [- o} end 
f(n,N) are given in Figs. 13 and 1}. 

Actually, over the interesting renge of parameters, the curve 
Py = (1 + out ) exp (-2u"] is not the best curve to use in (II.16) to give 
the best numerical results in approximating (II.14). It is possible to 
give this curve an average correction so as ta obtain better results over 
the interesting range of n and N. 

Such a corrected curve, to be used with the scale factor f(n,N) to 
approximate (II.14), is given in Fig. 13. The use of this curve with 
f(n,N) gives agreement with (II.14) to within about seven per cent in 

6< 12 


= for P. rol: 10 en 310 
Ro D 


better as Py gets larger, and is practically exact for Py - 50. 


» and 1 Sn $1000. the agreement gets 
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‘CASE | 
The greet density function f(v) is a polynomial multiplied by 
exp WA r 2). Except for small N, the polynomial is too long to be 
of much use for computation. Hence an Edgeworth series seems to be the 
best way to compute Py except for small N. | 
For N = 1, the exect formula is 


ahd 2 =p as 
Py = (a + 2) ls Bi hs ry | exp ee: (Ir.18) 


For N moderately large » the Edgeworth series is given by (II.11) with, 


letting B = 1 + 3 


> N(1 + x) 
ie 
N(2p~ - 1) 
3 
i oB~ -~ 1 
ett ayy (2p pl? 
4 
1 (2p" - 1) 
(e = 
4 IN (262 a 1° 
3 2 
1 2 1 
Cp = (II.19) 
6 BN (og? L 1y° 


R 
Curves of PD vB Ro for this case are given in Figs. 15 - 17. 


III, DERIVATION OF FORMULAS FOR P 


The characteristic function for the sum, v, of N normalized pulses 
of signal plus noise, with constant signal-to-noise power ratio x, for a 
square law detector, is* 


< a 
C(p) = co [-m (3) C3IE oT} 
(p+1) 


If x is now a random variable with density function w(x,x) » and x is 
constant for each group of N pulses but fluctuates independently from group 


to group, then the characteristic function for the sum of N pulses is, since 


w(x,x) =O for x40, _ Set g 
| C(p) = wii w(x,x) C(p) dx (III.2) 
0 
CASE 1 
- lL -x/x > 
Here w(x,x) = =e (x 2 0) 
therefore 
C(p) = neg (ILI.5) 
(l+p)"” [rep(2+ne)] . 
For N = 1, the density function ig(*) 
f(v) = s exp |r| (v 20) (III.4) 


For N >1, the density function can be found from pair 581.7 of Ref. 4%: 


uD v 
f(v) ie (2 + ar 5 Ns 9 ——-—_—, N-2 ool 4 (III. 5) 
Nx ab 1+Nx 
0 + ng) V-1 
therefore, for N =1, 


* In the notation used here, if f(v) is the density function and 
V(s) = fecrrent*” dv then p = is and C(p) = [(s). 
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P. <= exp ae (N = 1) (III.6) 


For N > 1, a simple way to get Py is: 


Gp) - eet = =p (1 + NE) G(p) 
+p 


e N-2_-v 
t(v) - Gaye = - (1+ Nx) f' (v) 


and therefore 


‘ t 
J ridarer | 2 ‘ r2| - (1 + Nx) £(¥,) 
0 N-1 


Since P, = 1 ah f£(v) dv, we have for N > 1: 
0 
N-1 z 
Pypei-I ear) ee Pie si Wis sy ee zh (TE.7) 
a Nx b+ =) N-1 1L + NX. 
Nxf | 
CASE 2 


The characteristic function for one pulse from a target fluctuating 


according to (I.1) is (from III.3 for N = 1): 


é = - 8 
 () mG Ane (III.8) 


Therefore the characteristic function for the sum of N pulses, 


fluctuating independently from pulse to pulse, is 


[e () | ; Bg pei REE TS (III.9) 
: =f. +p04 2)" 
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The density function is 


ely) « San aaa vith | exp ees (IIT.10) 
(14k) (N-1)! 
and therefore 
Py = i-t cha (III.11) 
(14%) YN 3 


Reference 3 contains tables of this function up to N - 1 = 50. Beyond 
(2) 


this, one can use an Edgeworth series to compute Pp: This can be done if 
one knows the cumulants K, for the density function f(v) (as well as the 
meanyise) 
In this case 
nai {a2 = 
K = N(-1) — [an(asp(1s))] 
dp px0 

or 


K, = N(n-1): (14x) (III.12) 


In addition, it is easily found that the mean of the sum of N pulses 
is N(l +x); the variance is (14x )*. Therefore the coefficients of the 


second, third, and fourth terms of the Edgeworth series are 


ee ee a 1 
Cx. a Bun? C), = iy? Ce " Ton (III.13) 
CASE 3 z 
1 hie ~2x/x 
Here w(x,x) = <5. (0 end therefore 
x 
= "i 
C(p) = (III.14) 


(14p )"-? [2 +p (2 + ey) 
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For N = 1 and N = 2, the density function is obtained from Ref. 4 as: 


Vv Vv 


Nel: f(v) c=] EER) ae 3 exp Per. 
Q+ 3] L+§ 1+3 

N= 2: f(v) = —+ eo | ty :| (III.15) 
(1 + x)* 


From Ref. 4, pair 581.1, after appropriate transformation, one finds 


that for N > 2, 


1 Nel -v Vv 
f(v) = 5 a aE v= 5 6 Fx | 2%, a =| (III.16) 
(N-1)! (2 ie x) Nx 


where Fi is the confluent hypergeometric function, * 


One may now use two identities concerning this hypergeometric function. ** 


a) F,(2,N,2) = (z +2-N) Fi (1,N,z)+N-1 


©) 4F,(1,N,z) =e” z Sa (y Seite eae he (112.17) 
| \N-i 
to get 
lL + Ned iy «y 
f(v) = 3 VoL . » N-2] exp iz 
(0. + 8) Q+% N-1 1+> 
N-1 
Cine rae | 
5 ae I Sy = » N-2/ exp| —Ss 
(. + 8) Q + ig) Va 1+> 
1 N- cE 
6g (III.18): 


* See Refs. 2 and 4. 
** Reference 2, mathematical appendix, pp 20-21. (In this reference there 
is a misprint in the printing of identity (b)). 
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If one assumes the false alarm probability to be negligible, and 


“1% 


to be >1, then for v >%, , the I factors are practically equal to one; 
also; the integral from ¥, to coof the third term in (IIT.18) can be 
neglected, for N>2. Therefore, integrating (III.18) from % to oo, 


one obtains for N > 2: 


N-2 2 
Fe (1 + &) 1+—he - ANAL) exp oe (111.19) 
+e Ll+> 


This is the exact formula for N = 1 and N = 2, as can be found directly 


from (III.15). 


CASE 4 


For one pulse (setting N = 1 in (III.14)) 


C, (2) = an SL, (III.20) 
ore) 


Therefore for the sum of N pulses, fluctuating independently from pulse to 


pulse, the characteristic function is 


(iD) |" = SOS as rh (11T.21) 
[é,(2)| ee ir IIT. 21 


The exact density function is obtainable and turns out to be a 


rol si 


polynomial mltiplied by exp [ = (ae x)| . Except for small N, this 
polynomial is too long to be useful for computational purposes. For 
fairly large N, one can use the Edgeworth series. For N =1, the exact 
formula is, of course, given by (III.15) with N= 1. 

The cumulants, obtained in the usual way, are 


x 


3 (III.22) 


K, = N(n-1)! [ 28" a 1] where B = 1 + 
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Also, mean = N(1l + x) 
variance = n(2p°=1) 
So,. the coefficients for the second, third, and fourth terms of the 


Edgeworth series are 


Pee es 28° «1 
5 osdl_ (28* 41)°/ 
a opt 1 
iN (28 21° 
Sas 
1 (26° -1) 
CC. = 
§ ~ IBN (a9? -1)9 
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(III.23) 


Appendix A 
DERIVATION OF EQUATION I1.8 


Let 
Nel 
fe 4 f 4 2 


Ou (2 ae 


then (referring to page 278 for definition of R, and R,) 


R, 4 R, 6 
rag aimee (=) ABE 4 
° ° 
80 
RAV kf nP 
2 1 D 
R er 2 
° 2 ft 
a D 
then, assuming <<l, 
fe 
R, to 
L 
where 
Nel 
penn 


E(n,N) = ny, 
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Case 3— approximate 
Pp=e-2u4 (1+2u4) 


—-—— — Corrected for better numerical results 
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